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We introduce the ring of Fermat reals, an extension of the real field containing 
nilpotent infinitesimals. The construction takes inspiration from Smooth Infinitesi- 
mal Analysis (SIA), but provides a powerful theory of actual infinitesimals without 
any need of a background in mathematical logic. In particular, on the contrary with 
respect to SIA, which admits models only in intuitionistic logic, the theory of Fermat 
reals is consistent with classical logic. We face the problem to decide if the product 
of powers of nilpotent infinitesimals is zero or not, the identity principle for poly- 
nomials, the definition and properties of the total order relation. The construction 
is highly constructive, and every Fermat real admits a clear and order preserving 
geometrical representation. Using nilpotent infinitesimals, every smooth functions 
becomes a polynomial because in Taylor's formulas the rest is now zero. Finally, we 
present several applications to informal classical calculations used in Physics: now 
all these calculations become rigorous and, at the same time, formally equal to the 
informal ones. In particular, an interesting rigorous deduction of the wave equation is 
given, that clarifies how to formalize the approximations tied with Hook's law using 
this language of nilpotent infinitesimals. 
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INTRODUCTION AND GENERAL PROBLEM 



Frequently in work by physicists it is possible to find informal calculations like 
1 v'^ 1 

1 + — - hu{x) = 1 - -/l44(x) (1) 



with explicit use of infinitesimals f/c <^ 1 or /i44(x) <^ 1 such that e.g. /i44(a;)^ = 0. For 



example Einstein 
equality sign ~) 



13l | wrote the formula (using the equality sign and not the approximate 



df 

f{x,t + T)=f{x,t)+T-^ix,t) (2) 



justifying it with the words "since t is very smalV; the formulas ^ are a particular case of 
the general ([2]) . Dirac jl^ wrote an analogous equality studying the Newtonian approxima- 
tion in general relativity. 

Using this type of infinitesimals we can write an equality^ in some infinitesimal neighbor- 
hood, between a smooth function and its tangent straight line, or, in other words, a Taylor's 
formula without remainder. 

There are obviously many possibilities to formalize this kind of intuitive reasonings, ob- 
taining a more or less good dialectic between informal and formal thinking, and indeed 
there are several theories of actual infinitesimals (from now on, for simplicity, we will say 
"infinitesimals" instead of "actual infinitesimals" as opposed to "potential infinitesimals"). 
Starting from these theories we can see that we can distinguish between two type of defi- 
nitions of infinitesimals: in the first one we have at least a ring R containing the real field 
M and infinitesimals are elements e E R such that — r < e < r for every positive standard 
real r G M>o- The second type of infinitesimal is defined using some algebraic property of 
nilpotency, i.e. = for some natural number n G N. For some ring R these definitions can 
coincide, but anyway they lead, of course, only to the trivial infinitesimal e = if i? = M. 

However these definitions of infinitesimals correspond to theories which are completely 
different in nature and underlying ideas. Indeed these theories can be seen in a more inter- 
esting way to belong to two different classes. In the first one we can put theories that need 
a certain amount of non trivial results of mathematical logic, whereas in the second one 
we have attempts to define sufficiently strong theories of infinitesimals without the use of 
non trivial results of mathematical logic. In the first class we have Non-Standard Analysis 
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(NSA) and Synthetic Differential Geometry (SDG, also called Smooth Infinitesimal Analy 



SIS, see e.g. Bell Kock H, Lavendhomme 



22|, Moerdijk and Reyes [2^), in the second 



2 lip. Levi-Civita fields (see Berz 



one we have, e.g., Weil functors (see Kriegl and Michor 

Shamseddine 25|), surreal numbers (see Conway j^, Ehresmann geometries over 



rings containing infinitesimals (see Bertram [6!|). More precisely we can say that to work in 
NSA and SDG one needs a formal control deeply stronger than the one used in "standard 
mathematics". Indeed to use NSA one has to be able to formally write the sentences one 
needs to use the transfer theorem. Whereas SDG does not admit models in classical logic, 
but in intuitionistic logic only, and hence we have to be sure that in our proofs there is no 
use of the law of the excluded middle, or e.g. of the classical part of De Morgan's law or of 
some form of the axiom of choice or of the implication of double negation toward affirmation 
and any other logical principle which is not valid in intuitionistic logic. Physicists, engineers, 
but also the greatest part of mathematicians are not used to have this strong formal control 
in their work, and it is for this reason that there are attempts to present both NSA and 
SDG reducing as much as possible the necessary formal control, even if at some level this 

n riR 

is technically impossible (see e.g. Henson [19], and Benci and Di Nasso [J, 15|| for NSA; Bell 



and Lavendhomme 



22| for SDG, where using an axiomatic approach the authors try to 



postpone the very difficult construction of an intuitionistic model of a whole set theory using 
Topos). 

On the other hand NSA is essentially the only theory of infinitesimals with a discrete 
diffusion and a sufficiently great community of working mathematicians and published results 
in several areas of mathematics and its applications, see e.g. Albeverio et al. SDG is 
the only theory of infinitesimals with non trivial, new and published results in differential 
geometry concerning infinite dimensional spaces like the space of all the diffeomorphisms of 
a generic (e.g. non compact) smooth manifold. In NSA we have only few results concerning 
differential geometry. Other theories of infinitesimals have not, at least up to now, the same 
formal strength of NSA or SDG or the same potentiality to be applied in several different 
areas of mathematics. 

Our main aim, of which the present work represents a first step, is to find a theory 
of infinitesimals within "standard mathematics" (in the precise sense explained above of a 
formal control more "standard" and not so strong as the one needed e.g. in NSA or SDG) 
with results comparable with those of SDG, without forcing the reader to learn a strong 
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formal control of the mathematics he/she is doing. Because it has to be considered inside 
"standard mathematics", our theory of infinitesimals must be compatible with classical logic. 

Concretely, the idea of the present work is to by-pass the impossibility theorem about 
the incompatibility of SDG with classical logic that forces SDG to find models within intu- 
itionistic logic. 

Another point of view about present theories of infinitesimals is that, in spite of the 
fact that frequently they are presented using opposed motivations, they lacks the intuitive 
interpretation of what the powerful formalism permits to do. For some concrete example 
in this direction, see Giordano jl6|. Another aim of the present work is to construct a new 
theory of infinitesimals preserving always a very good dialectic between formal properties 
and intuitive interpretation. 

More technically we want to show that it is possible to extend the real field adding 
nilpotent infinitesimals, arriving at an enlarged real line 'M, by means of a very simple 
construction completely inside "standard mathematics". Indeed to define the extension 'M D 
M we shall use elementary analysis only. To avoid misunderstandings is it important to clarify 
that the purpose of the present work is not to give an alternative foundation of differential 
and integral calculus (like NSA), but to obtain a theory of nilpotent infinitesimals as a first 
step for the foundation of a smooth (C°°) differential geometry. For some preliminary results 
in this direction, see Giordano [l 



II. MOTIVATIONS FOR THE NAME "FERMAT REALS" 

It is well known that historically two possible reductionist constructions of the real field 
starting from the rationals have been made. The first one is Dedekind's order completion 
using sections of rationals, the second one is Cauchy's metric space completion. Of course 
there are no historical reason to attribute our extension 'M D M of the real field, to be 
described below, to Fermat, but there are strong motivations to say that, probably, he 
would have liked the underlying spirit and some properties of our theory. For example: 

1. a formalization of Fermat's infinitesimal method to derive functions is provable in our 
theory. We recall that Fermat's idea was, roughly speaking and not on the basis of an 
accurate historical analysis which goes beyond the scope of the present work (see e.g. 



Edwards [U|, Eves [l^), to suppose first h ^ 0, to construct the incremental ratio 

fix + h)- fix) 
h 

and, after suitable simplifications (sometimes using infinitesimal properties), to take 
in the final result h = 0. 

2. Fermat's method to find the maximum or minimum of a given function fix) at x = a 
was to take e to be extremely small so that the value of /(x + h) was approximately 
equal to that of /(x). In modern, algebraic language, it can be said that fix+h) = fix) 
only if h"^ = 0, that is if e is a first order infinitesimal. Fermat was aware that this 
is not a "true" equality but some kind of approximation (ibidem). We will follow a 
similar idea to define 'M introducing a suitable equivalence relation to represent this 
equality. 

3. Fermat has been described by Bell ^ as "the king of amateurs" of mathematics, and 
hence we can suppose that in its mathematical work the informal/intuitive part was 
stronger with respect to the formal one. For this reason we can think that he would 
have liked our idea to obtain a theory of infinitesimals preserving always the intuitive 
meaning and without forcing the working mathematician to be too much formal. 

For these reason we chose the name "Fermat reals" for our ring 'R (note: without the 
possessive case, to underline that we are not attributing our construction of 'M to Fermat). 



III. DEFINITION AND ALGEBRAIC PROPERTIES OF FERMAT REALS: THE 

BASIC IDEA 

We start from the idea that a smooth (C°°) function / : 'M — > 'M. is actually equal to 
its tangent straight line in the first order neighborhood e.g. of the point x = 0, that is 

yheD: fih) = fiO) + h-fiO) (3) 

where D is the subset of 'M which defines the above-mentioned neighborhood of x = 0. 
The equality ([3]) can be seen as a first-order Taylor's formula without remainder because 
intuitively we think that h'^ = for any h G D (indeed the property h'^ = defines the 
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first order neighborhood of x = in *R). These almost trivial considerations lead us to 
understand many things: 'M must necessarily be a ring and not a field because in a field 
the equation = implies h = 0; moreover we will surely have some limitation in the 
extension of some function from R to *M, e.g. the square root, because using this function 
with the usual properties, once again the equation = implies \h\ = 0. On the other 
hand, we are also led to ask whether (JS]) uniquely determines the derivative /'(O): because, 
even if it is true that we cannot simplify by h, we know that the polynomial coefficients of 
a Taylor's formula are unique in classical analysis. In fact we will prove that 

3! m e M V/i e D : f{h) = /(O) + h ■ m (4) 

that is the slope of the tangent is uniquely determined in case it is an ordinary real number. 
We will call formulas like dH) derivation formulas. 

If we try to construct a model for (J!]) a natural idea is to think our new numbers in *M. 
as equivalence classes [h] of usual functions h : M — > M. In this way we may hope both 
to include the real field using classes generated by constant functions, and that the class 
generated by h(t) = t could be a first order infinitesimal number. To understand how to 
define this equivalence relation we have to think at (I3|) in the following sense: 

/(Mt))~/(0) + Mt)-/'(0), (5) 

where the idea is that we are going to define ~. If we think h{t) "sufficiently similar to f , 
we can define ~ so that ((51) is equivalent to 

^.^ f{h{t)) - /(O) - h{t) ■ f (0) ^ Q 

i-*0+ t ' 

that is 

X ~ y : <^==^ lim — — ^ = 0. (6) 

In this way ^ is very near to the definition of differentiability for / at 0. 

It is important to note that, because of de L'Hopital's theorem we have the isomorphism 

C^(R,R)/~ - RH/(x), 

the left hand side is (isomorphic to) the usual tangent bundle of R and thus we obtain 
nothing new. It is not easy to understand what set of functions we have to choose for x, y 
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in ([6l) so as to obtain a non trivial structure. The first idea is to take continuous functions 
at t = 0, instead of more regular ones like C^-functions, so that e.g. hk{t) = \t\^^'' becomes a 
k-th order nilpotent infinitesimal ~ 0); indeed for almost all the results presented in 

this article, continuous functions at t = work well. However, only in proving the non-trivial 
property 

(Vx G 'M : X ■ fix) =0) =^ Vx G 'R : /(x) = (7) 

we can see that it does not suffice to take continuous functions at t = 0. To prove ([7]) the 
following functions turned out to be very useful: 

Definition 1. If x : ]R>o — > M, then we say that x is nilpotent iff \x(t) — x(0)|'^ = o{t) as 
t ^ 0^, for some k eN. Af will denote the set of all the nilpotent functions. 

E.g. any Hoelder function |x(t) — x(s)| < c ■ |t — s|" (for some constant a > 0) is nilpotent. 
The choice of nilpotent functions instead of more regular ones establish a great difference of 
our approach with respect to the classical definition of jets (see e.g. Brocker [8|, Golubitsky 



and Guillemin [l7||), that ^ may recall. 

Another problem necessarily connected with the basic idea ([3]) is that the use of nilpotent 
infinitesimals very frequently leads to consider terms like h^^-. ■ --h^n ■ For this type of products 
the first problem is to know whether ■ . . . ■ h"^ 7^ and what is the order k of this new 
infinitesimal, that is for what k we have [h!^ ■ . . . ■ k"^)^ ^ but {W^ ■ . . . ■ = 0. We 

will have a good frame if we will be able to solve these problems starting from the order of 
each infinitesimal hj and from the values of the powers ij G N. On the other hand almost all 
the examples of nilpotent infinitesimals are of the form h{t) = t", with < a < 1, and their 
sums; these functions have also great properties in the treatment of products of powers. It 
is for these reasons that we shall focus our attention on the following family of functions 
X : ]R>o — > M in the definition ([6]) of ~: 

Definition 2. We say that x is a little-oh polynomial^ and we write x G Mo[i] iff 

1. X : M>o — > M 

2. We can write 

k 

x^ = r + ^ai ■ t"' + o{t) as t ^ 0+ 
1=1 
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for suitable 

ken 

r,ai, . . . ,ak eR 
ai, . . . ,ak G M>o 

Hence a little-oh polynomial x e Mo[^] is a polynomial function with real coefficients, in 
the real variable t > 0, with generic positive powers of t, and up to a little-oh function as 

Remark 3. In the following, writing Xt = yt + o{t) as t — * 0+ we will always mean 

lim — — ^ = and xq = uq. 

t^o+ t 

In other words, every little-oh function we will consider is continuous as t — >^ 0+. 
Example. Simple examples of little-oh polynomials are the following: 

1. Xt = 1 + t + + + o{t) 

2. Xt = r Wt. Note that in this example we can take k 
void sequence of reals, that is the function a = a : 
X of reals as a function x : {1, . . . ,n} — > M. 

3. Xt = r + o{t) 

IV. FIRST PROPERTIES OF LITTLE-OH POLYNOMIALS 

Little-oh polynomials are nilpotent: 

First properties of little-oh polynomials are the following: if = r + ^*L^ aj-t"'+Oi(t) as 
t ^ 0+ and yt = s + Ef=i /^i "^'^ +02(t), then {x+y) = r + s + ^ti «i-^"" + Ef=i /^r^'^ +03(t) 
and {x ■y)t = rs + Ya=i -^^i " + SjLi ■ t^' + Ya=i Y.f=i CdPj ■ ^""'t^' + Oi{t), hence the 
set of little-oh polynomials is closed with respect to pointwise sum and product. Moreover 
little-oh polynomials are nilpotent (see Definition [T|) functions; to prove this we firstly prove 
that the set of nilpotent functions A/" is a subalgebra of the algebra M"^ of real valued 
functions. Indeed, let x and y be two nilpotent functions such that \x — a;(0)|'^ = Oi(t) and 



0, and hence a and a are the 
M, if we think of an n-tuple 
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\y — y{0)\^ — 02{t), then we can write x-y — x{0)-y{0) — x -[y — y{0)]+y{0) - [x — x{0)], so that 
we can consider \x ■ [y - y(0)]|'= = {xl'' ■ \y - y(0)|*^ = 1x1*= • oi{t) and M^|i(*) ^ q as i ^ 0+ 
because \x\'' |x(0)|'^, hence x ■ [y — y{0)] G Af. Analogously y{0) ■ [x — x{0)] G Af and 
hence the closure of A/" with respect to the product follows from the closure with respect to 
the sum. The case of the sum follows from the following equalities (where we use xt := x{t), 
u :— X — xo, V :— y — yo, \ut\'' — oi{t) and \vt\^ — 02{t) and we have supposed k > N): 

= Oi(t), v'' = 02 (t) 
k 



Vi = 0, . . . , A; : 



i=0 ^ ^ 



t tk ■ t k \ ^ J \ ^ 

Xow we can prove that Mo[t] is a subalgebra of A/". Indeed every constant r G M and every 
power t""^ are elements of Af and hence r + Yli=i ' ^ it remains to prove that if 

y & Af and w — o{t), then y+w G Af, but this is a consequence of the fact that every little-oh 
function is trivially nilpotent, and hence it follows from the closure of jV' with respect to the 
sum. 



Closure of little-oh polynomials with respect to smooth functions: 

Now we want to prove that little-oh polynomials are preserved by smooth functions, that 
is if a; G MqM and / : M — > M is smooth, then / o x G Mo[t]. Let us fix some notations: 

k 

= r + ^ ctj • t"' + w{t) with w{t) = o{t) 

i=l 

h{t) := x{t) - x{0) Vt G M>o 

hence Xt — x{0) + ht — r + ht- The function t i— > h{t) — Yli=i ' + '^{^) belongs to 
^o\t\ C jV so we can write \h\^ — o{t) for some A^" G N and as i — > 0+. From Taylor's 
formula we have 

f{xt) = /(r + = /(r) + ■ K + fi^t) = /(r + h^) (8) 



1=1 



^ t{'i')( \ 

M + E ■ + o{h^) (9) 



i=l 
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But 

\o{hf)\ _ |o(/.f)| \hll 

\t\ K\ ' \t\ ^ 

hence o{hf) = o{t) G Mo[^]- From this, the formula (l8|), the fact that h G M.o[t] and using the 
closure of little-oh polynomials with respect to ring operations, the conclusion f o x E M.o[t] 
follows. 



V. EQUALITY AND DECOMPOSITION OF FERMAT REALS 

Definition 4. Let x, y E Mo[^], then we say that a; ~ ?/ or that a; = ?/ in 'M iff x{t) = 
y{t) + o{t) as t 0'^ . Because it is easy to prove that ~ is an equivalence relation, we can 
define 'M := MoM/ ~, i-e. *M is the quotient set of M.o[t] with respect to the equivalence 
relation ~. 

The equivalence relation ~ is a congruence with respect to pointwise operations, hence 'R 
is a commutative ring. Where it will be useful to simplify notations we will write "x = y 
in "M" instead of x ~ ?/, and we will talk directly about the elements of M.o[t] instead of 
their equivalence classes; for example we can say that x = y in 'M and z = w in *M. imply 
x + z = y + w in "M. 

The immersion of M in *M is r i — > f defined by f(t) := r, and in the sequel we will always 
identify M with M, which is hence a subring of *]R. Conversely if x G "M then the map 
°(— ) : X E *M. \—>- °x = x{0) G M, which evaluates each extended real in 0, is well defined. 
We shall call °(— ) the standard part map. Let us also note that, as a vector space over the 
field M we have dim^ 'M = oo, and this underlines even more the difference of our approach 
with respect to the classical definition of jets. Our idea is instead more near to NSA, where 
standard sets can be extended adding new infinitesimal points, and this is not the point of 
view of jet theory. 

With the following theorem we will introduce the decomposition of a Fermat real x G *M, 
that is a unique notation for its standard part and all its infinitesimal parts. 

Theorem 5. If x E *M, then there exist one and only one sequence 

{k,r,ai, ...,ak,ai,.. .,ak) 

such that 

kEN 
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r, ai, . . . ,afc,ai, . . . ,afc G M 

and 

k 

L x = r + Y,ai- in 'M 

1=1 

2. < ai < a2 < ■ ■ ■ < Ofc < 1 
5. ^ Vi = 1, . . . , 

In this statement we have also to include the void case /c = and a = a : — > M. Obviously, 
as usual, we use the definition Y^i=i h = for the sum of an empty set of numbers. As we 
shall see, this is the case where x is a standard real, i.e. x G M. 

In the following we will use the notations t° := dti/a := [t G M>o i-^ G M]^ G 'M so that 
e.g. dt2 = t^^'^ is a second order infinitesimal. In general, as we will see from the definition 
of order of a generic infinitesimal, dta is an infinitesimal of order a. In other words these 
two notations for the same object permit to emphasize the difference between an actual 
infinitesimal dta and a potential infinitesimal t^/": an actual infinitesimal of order a > 1 
corresponds to a potential infinitesimal of order ^ < 1 (with respect to the classical notion 
of order of an infinitesimal function from calculus, see e.g. Prodi Silov ^^). 

Remark 6. Let us note that dta ■ dt^ = dt_ab_, moreover dt" := (dta)"' = dts. for every 

a + b " ^ ' O! 

a > 1 and finally dta = for every a < 1. E.g. dtL"'^^ = for every a G M>o, where [a] G N 
is the integer part of a, i.e. [a] < a < [a] + 1. 

Existence proof: 

Since x G ]Ro[t], we can write Xt = r + Yl'i=i '^j ' + *^(^) as t ^ 0+, where r, ctj G M, 
Oj G M>o and /c G N. Hence x = r + Yl^=i ' '^^^ purpose is to pass from 

this representation of x to another one that satisfies conditions [H [2] and [3] of the statement. 
Since if Oj > 1 then ■ t"' = in 'M, we can suppose that ai < 1 for every i = 1, . . . ,k. 
Moreover we can also suppose > for every i, because otherwise, if = 0, we can replace 
r G M by r + I = 0, ^ = 1, • • • , k}. 

Now we sum all the terms t"'^ having the same Oj, that is we can consider 

di := ^{oj \ aj = tti , j = 1, . . . , k} 

so that in 'M. we have 

X = r + ^di-t"-' 
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where / C {1, . . . , fc}, {oj \ i E 1} = {a, . . . , a^} and ^ aj for any i, j E I with i ^ j. 
Neglecting a, if a, = and renaming ai, for i G /, in such a way that < aj if i, j G / with 
? < j, we obtain the existence result. Note that if a; = r G M, in the final step of this proof 
we have / = 0. 

Uniqueness proof: 

Let us suppose that in 'M we have 

k N 

X = r + '^ai-t^' = s + ^(3j ■ t^^ (10) 

i=i j=i 

where ai, j3j, ai and hj verify the conditions of the statement. First of all °x = x{0) = r = s 
because ai, bj > 0. Hence ait"^ — f3it^^ + Zli '^j ' ~ Z^j ' = o{t). By reduction to the 
absurd, if we had ai < 6i, then collecting the term f^^ we would have 

* 3 

In ( ITTi l we have that I3it^^~°-^ — > for t — > 0+ because oi < hi by hypothesis; oii-t°'^~"'^ 
because Oi < Oj for i = 2, . . . , A;; /3j ■ — > because ai < 6i < bj for j = 2, . . . , A^, 
and finally )f:^~'*i is limited because ai < 1. Hence for t — 0+ we obtain ai = 0, which 
conflicts with condition [3] of the statement. We can argue in a corresponding way if we had 
bi < ai. In this way we see that we must have Oi = bi. From this and from equation ([TTIl 
we obtain 

i j 

and hence for t — > 0^ we obtain ai = Pi. We can now restart from (fT2ll to prove, in the 
same way, that 02 = &2, ^2 = /?2, etc. At the end we must have k = N because, otherwise, 
if we had e.g. < A^, at the end of the previous recursive process, we would have 

TV 

j=k+i 

From this, collecting the terms containing t^'^+i, we obtain 

1^"^+'-^ ■ [(3k+i + pk+2 ■ 1^"+^-^"+^ + . . . + . t/5iv-A+i] ^ 0. (13) 

In this sum (3k+j " as t ^ 0"*", because 6fc+i < 6^+^ for j > 1 and hence 

(3k+i + (3k+2 ■ + . . . + . tf^N-Pu+i ^ ^^^^ -L 0, so from (HI) we get t^^+^-^ 0, 

that is bk+i > 1, in contradiction with the uniqueness hypothesis bk+i < 1. 
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Let us note explicitly that the uniqueness proof permits also to affirm that the decom- 
position is well defined in 'M, i.e. that if x = y in 'M, then the decomposition of x and the 
decomposition of y are equal. 

■ 

On the basis of this theorem we introduce two notations: the first one emphasizing the 
potential nature of an infinitesimal x G 'R, and the second one emphasizing its actual 
nature. 

Definition 7. If x G 'M, we say that 

k 

X = r + ai ■ t"' is the potential decomposition (of x) (14) 

i=l 

iff conditions [H, El , and[3l of Theorem [5] are verified. Of course it is implicit that the symbol 
of equality in ffT4ll has to be understood in 'R. 

For example a; = 1 + t^^^ + t^/^ + t is a decomposition because we have increasing powers 
of t. The only decomposition of a standard real r G R is the void one, i.e. that with k = 
and a = a : — > R; indeed to see that this is the case, it suffices to go along the existence 
proof again with this case a; = r G R (or to prove it directly, e.g. by contradiction). 

Definition 8. Considering that = dti/a, we can also use the following notation, empha- 
sizing more the fact that a; G *R is an actual infinitesimal: 

k 

x = °x + J2°^i- (15) 

j=i 

where we have used the notation °Xi := and 6j := l/ctj, so that the condition that uniquely 
identifies all 6j is 6i > 62 > ■ ■ ■ > bk > 1. We call (fTSll the actual decomposition of x or 
simply the decomposition of x. We will also use the notation d*x := °Xi ■ dt{,. (and simply 
da; := d^a;) and we will call °Xi the i-th standard part of x and d*x the i-th infinitesimal 
part of X or the i-th differential of x. So let us note that we can also write 

a; = °a; + d*a; 

i 

and in this notation all the addenda are uniquely determined (the number of them too). 
Finally, if A; > 1 that is if a; G *R \ R, we set uj{x) := hi and uji{x) := hi. The real 
number uj{x) = hi is the greatest order in the actual decomposition (fTSj) . corresponding to 
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the smallest in the potential decomposition (fT4l) . and is called the order of the Fermat real 
X G *M. The number uji{x) = hi is called the i-th order of x. If a; G M we set uj{x) := and 
d*x := 0. Observe that in general uj{x) = ^^(da;), d(dx) = da; and that, using the notations 
of the potential decomposition ((71), we have uj{x) = 1/ai. 

Example. If x = 1 + t^^^ + t^/^ + t, then °x = 1, dx = dt^ and hence x is a third order 
infinitesimal, i.e. u{x) = 3, d^x = dt2 and d^x = dt; finally all the standard parts are 
°Xi = 1. 

VI. THE IDEALS Dk 

In this section we will introduce the sets of nilpotent infinitesimals corresponding to 
a k-th order neighborhood of 0. Every smooth function restricted to this neighborhood 
becomes a polynomial of order k, obviously given by its k-th order Taylor's formula (without 
remainder). We start with a theorem characterizing infinitesimals of order less than k. 

Theorem 9. //x G 'M and k G N>i, then x'^ = m 'M if and only if°x = and uj{x) < k. 

Proof: If x'^ = 0, then taking the standard part map of both sides, we have "(x'^) = (°x)^ = 
and hence °x = 0. Moreover x'' = means x^ = o{t) and hence {^rfw)^ ~^ and -jjfjr — > 0. 
We rewrite this condition using the potential decomposition x = Yl'i=i ' ^""^ of ^ (note that 
in this way we have u;(x) = ^) obtaining 

lim V (Xi ■ f'-T^ = = lim t'^^"^ ■ \ai + ■ t"^'"' + • • • + ■ t""-'''] 

i 

But ai + ^2 ■ r^-ai _^ . . . _^ Q,^ . fcik-ai _^ q,^ _^ hcucc wc must have that t"^^^ 0, and 
so oi > |, that is uj{x) < k. 

Vice versa if °x = and lj{x) < fc, then x = Yli=i ' + '^(^)i 

hm -TT = hm > ■ + lim — ■ t^'^ 

i 

But t^'T^ because k> 1 and t°-^~T^ 0+ because — < — = u!(x) < k and hence x'^ = 
in 'R. 

■ 

If we want that in a k-th order infinitesimal neighborhood a smooth function is equal to its 
k-th Taylor's formula, we need to take infinitesimals which are able to delete the remainder. 
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that is, such that = 0. The previous theorem permits to extend the definition of the 
ideal Dk to real number subscripts instead of natural numbers k only. 

Definition 10. If a G M>o U {oo}, then 

:= {a; G 'M I °x = 0, uj{x) < a + 1} 
Moreover, we will simply denote Di hy D. 

1. If X = dts, then uj{x) = 3 and x G -D3. More in general dtk G Da if and only if 
Lu{ dtk) = k < a + 1. E.g. dtk e D if and only if 1 < A; < 2. 

2. Doo = [jaDa = {x e'R\°x = 0} is the set of all the infinitesimals of 'M. 

3. Dq = {0} because the only infinitesimal having order strictly less than 1 is, by defini- 
tion of order, x = (see the Definition [8j). 

The following theorem gathers several expected properties of the sets Da and of the order 
of an infinitesimal uj{x). 

Theorem 11. Let a, b E M>o and x, y E Doc, then 

1. a<h =^ Da^Db 

2. X e D^(^) 

5. a G N =^ Da = {x e'R\ x''^^ = 0} 

4. X e Da =^ X^''^^^ = 

5. X G -Doo \ {0} and k = [uj{x)] =^ x E Dk \ -D/t-i 

6. d{x ■ y) = dx ■ dy 

1 1 1 

7. x-y^O =^ — = — - + — - 

uj{x-y) uj{x) uj[y) 

8. x + y^O =^ Lj{x + y) = uj{x) y Lj{y) 

9. Da is an ideal 
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In this statement if r G M, then \r] is the ceiling of the real r, i.e. the unique integer \r] E Z 
such that \r] — 1 < r < \r]. Moreover if r, s G M, then r V s := max(r, s). 

Property [7} of this theorem cannot be proved substituting the ceiling \a\ with the 
integer part [a]. In fact if a = 1.2 and x = dt2.i, then uj{x) = 2.1 and [a] + 1 = 2 so that 

^[a]+i = 3,2 ^ 7^ in 'M, whereas [a] + 1 = 3 and = dt^ = 0. 

2 ' ' 3 

Finally let us note the increasing sequence of ideals/neighborhoods of zero: 

{0} = DoCD = DiCD2C---CDkC---CD^. (16) 

Because of ( |T6l) and of the property dta = if a < 1, we can say that dt is the smallest 
infinitesimals and dt2, dts, etc. are greater infinitesimals; as we will see, this agree to 
corresponding order properties of these infinitesimals. 



VII. PRODUCTS OF POWERS OF NILPOTENT INFINITESIMALS 



In this section we will introduce some instruments that will be very useful to decide 
whether a product of the form h]^ ■ . . . ■ h^^ , with hk G -Doo \ {0}, is zero or whether it 
belongs to some D^. Generally speaking this problem is not trivial in a ring (e.g. in SDG 



there is not an effective procedure to decide this problem, see e.g. Lavendhomme 
its solutions will be very useful in the proofs of infinitesimal Taylor's formulas. 

Theorem 12. Let /ii, . . . , /i„ G -Doo \ {0} and zi, . . . , G N, then 

n 

1 " 



22|) and 



Proof: Let 



r=l 

be the potential decomposition of hk for k = 1, . . . ,n. Then by Definition [7] of potential 
decomposition and Definition [8] of order, we have < a^i < ak2 < ■ ■ < auN^ ^ 1 ^ind 
ik ■= uj{hk) = hence < aur for every r = 1, . . . .N^. Therefore from (fTTll . collecting 
the terms containing we have 

hk = t'/^^ ■ {aki + afcar"'"'/"'" + ■ ■ ■ + akN.t'"'^^-^^^^) 
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and hence 



'i'l \ I in / 1 

/ 1 1 \ «n 

• . . ■ (ani + an2f"'~^ + ■■■ + anNj"'-''"'^) (18) 

Hence if 22k ^ "'" have that t^i J" = in 'M, so also ■ . . . ■ = 0. Vice versa if 
/i*^^ ■ . . . ■ /i^" = 0, then the right hand side of (fTSl ) is a o(t) as t ^ 0+, that is 

lL+...+in_i / ,ai2-^ , , .aiN,--y^ 

Jn . ( Q,^^ _^ Q;^2t ^1 H haiATjt ^ ^1 ) 

/ 1 1 \ 

. . . • (a„l + an2t""'"^ + ■ ■ ■ + «nAr„t'"'^""^ j ^ 
/ 1 1 \ ^fc ■ 

But each term (aui + «fc2^"''^ + ■ ■ ■ + dkHk^'^'''^^ ^'^ ) ^ 7^ so, necessarily, we must 

have ^ + - + — 1>0, and this concludes the proof of d 

To prove O it suffices to apply recursively property [3 of Theorem [TTl 



Example 13. cu(dC, ■ • • • ■ dCj-^ = = t d^, • . . . ■ d^ = if and 

only if Ylik ^ ^ ^-S- dt ■ /i = for every /i G -Doo- 



The following corollary gives a necessary and sufficient condition to have ■ . . . ■ h^^ G 
D,\{^}. 

Corollary 14. In the hypotheses of the previous Theorem[TM let p G M>o, then we have 

^ n 

h{^-...-h:-eD,\{0} ^ — T<E^^l 

V+l ^ cuihk) 



Let h, k E D; because in this case J2k jf+T ~ i + I ~ -'- always have 

h-k = 0. (19) 

This is a great conceptual difference between Fermat reals and the ring of SDG, where, 
not necessarily, the product of two first order infinitesimal is zero. The consequences of 
this property of Fermat reals arrive very deeply in the development of the theory of Fermat 
reals, forcing us, e.g., to develop several new concepts if we want to generalize the derivation 
brmula (111) to functions defined on infinitesimal domains, like / : D — > 'M. (see Giordano 
IgI). We only mention here that looking at the simple Definition [U the equality (fTOll has 
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an intuitively clear meaning, and it is to preserve this intuition that we keep this equality 
instead of changing completely the theory toward a less intuitive one. 
Let us note explicitly that the possibility to prove these results about products of powers 
of nilpotent infinitesimals is essentially tied with the choice of little-oh polynomials in the 



definition of the equivalence relation ~ in Definition O Equally effective and useful resu 
are not provable for the more general family of nilpotent functions (see e.g. Giordano 
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VIII. IDENTITY PRINCIPLE FOR POLYNOMIALS AND INVERTIBLE 

FERMAT REALS 

In this section we want to prove that if a polynomial oq + ciix + 02^^ + ■ ■ ■ + a„x" of 'M 
is identically zero, then = for all /c = 0, . . . , n. To prove this conclusion, it suffices to 
mean "identically zero" as "equal to zero for every x belonging to the extension of an open 
subset of M". Therefore we firstly define what this extension is. 

Definition 15. If U is an open subset of M", then 'U := {x G 'M" | °x E U}. Here with the 
symbol 'R" we mean := 'R x . . . . x'M. 

The identity principle for polynomials can now be stated in the following way and proved 
in standard manner using Vandermonde matrices. 

Theorem 16. Let Oq, . . . , a„ G 'M and U be an open neighborhood of m R such that 

ao + aix + a2X^ H h a„x" = in'R Vs G 'U (20) 

Then 

ao = cti = ■ ■ ■ = a„ = in 'R 

Now, we want to see more formally that to prove ([3]) we cannot embed the reals M into a 
field but only into a ring, necessarily containing nilpotent element. In fact, applying ([3]) to 
the function f{h) = for h E D, where C 'R is a given subset of 'R, we have 

f{h) = h^ = /(o) + h ■ /'(o) = WheD. 

Where we have supposed the preservation of the equality /'(O) = from R to "R. In other 
words, if D and f{h) = verify ((31), then necessarily each element h E D must be a new 
type of number whose square is zero. 
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Because we cannot have property ([3]) and a field at the same time, we need a sufficiently 
good family of cancellation laws as substitutes. The simplest one of them is also useful to 
prove the uniqueness of 

Theorem 17. If x ^ 'M is a Fermat real and r, s G M are standard real numbers, then 

X ■ r = X ■ s m *]R and x 7^ =^ r = s 
Proof: From the Definition [4] of equality in 'M and from x ■ r = x ■ s we have 

hm = 0. 

t^o+ t 

But if we had r 7^ s this would implies limt^o+ = 0? that is a; = in 'M and this contradicts 
the hypothesis x ^ 0. 

m 

The last result of this section takes its ideas from similar situations of formal power series 
and gives also a formula to compute the inverse of an invertible Fermat real. 

Theorem 18. Let x = °x + Y17=i ' ^^ai be the decomposition of a Fermat real x E *M. 
Then x is invertible if and only if °x 7^ 0, and in this case 

j=0 \ i=l / 

In the formula (i2T]l we have to note that the series is actually a finite sum because any dta. 
is nilpotent, e.g. (1 + dt2)"^ = 1 — dt2 + dt^ — dt^ + ■ ■ ■ = 1 — dt2 + dt because dtl = 0. 

Proof: U X ■ y = 1 for some y G 'M, then, taking the standard parts of each side we 
have °x ■ °y = 1 and hence °x 7^ 0. Vice versa let y := °x~^ ■ X]j=o(~l)"' ' (X^i ^ dta^Y and 

h := X — °x = J2i dta.^ G -Dqo so that we can also write 

+00 ij 

j=0 

But /i G "M is a little-oh polynomial with h{0) = 0, so it is also continuous, hence for a 
sufficiently small 5 > we have 

ytei-5,6): <1. 

Therefore 

Vt G i-S, 5): yt = — .il + ^) = = - 

^ ^ °x \ °xj °x + ht Xt 

From this equality and from Definition [4] it follows x ■ y = \ in *]R. 
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IX. THE DERIVATION FORMULA 

In this section we want to give a proof of dH) because it has been the principal motivation 
for the construction of the ring of Fermat reals "K. Anyhow, before considering the proof of 
the derivation formula, we have to extend a given smooth function / : M — > M to a certain 
function 'f -.'R — > 'M. 

Definition 19. Let A be an open subset of M", / : A — > M a smooth function and x ^ 'A 
then we define 

'fix) ■■= f ox. 

This definition is correct because we have seen that little-oh polynomials are preserved by 
smooth functions, and because the function / is locally Lipschitz, so 



f{xt) - fiVt) 


<K- 


xt - Vt 


Vt e 


t 




t 





for a sufficiently small S and some constant K, and hence if x = y in 'R, then also ' f{x) = 
'f{y) in 'M. 

The function '/ is an extension of /, that is 

•/(r) = /(r) in 'M Vr G M, 

as it follows directly from the definition of equality in 'M (i.e. Definition SI), thus we can still 
use the symbol f{x) both for x G 'R and x G M without confusion. After the introduction of 
the extension of smooth functions, we can also state the following useful elementary transfer 
theorem for equalities, whose proof follows directly from the previous definitions: 

Theorem 20. Let A be an open subset of W, and r, a : A — > M be smooth functions. 
Then it results 

Vx G M : •r(x) = V(x) 

Vr G A : r(r) = <j{r). 
Now we will prove the derivation formula (l4|). 

Theorem 21. Let A be an open set m R, x G A and f : A — > R a smooth function, then 

3! m G R MheD : f{x + h) = /(x) + h ■ m. (22) 
In this case we have m = f'{x), where f'{x) is the usual derivative of f at x. 
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Proof: Uniqueness follows from the previous cancellation law Theorem [T7l indeed if mi G M 
and 7712 G both verify (l22l) . then h ■ rrii = h ■ 1712 ior every h E D. But there exists a non 
zero first order infinitesimal, e.g. dt E D, so from Theorem f flTl ) it follows mi = m2. 

To prove the existence part, take h E D, so that /i^ = in 'M, i.e. /i^ = o{t) for t ^ 0+. 
But / is smooth, hence from its second order Taylor's formula we have 

fix + h) = fix) + h, . fix) + I ■ fix) + oihf) 



But 



so 



oih^) oih^) 

= ■ ^ ^ for t 

t hj t 



^ . f" ix) + oih]) = o,it) fort^O^ 



and we can write 



fix + ht) = fix) + ht ■ fix) + Oi(t) for t ^ 0+ 

that is 

fix + h) = fix) + h ■ fix) in 'R 
and this proves the existence part because fix) E M. 

■ 

For example = 1 + h, sin(/i) = h and cos(/i) = 1 for every h E D. 
Analogously we can prove the following infinitesimal Taylor's formula. 

Lemma 22. Let A be an open set in R'^, x E A, n E N>o and f : A — > R a smooth 
function, then 

^ hi f)\j\f 
yhEDt: f{^ + h)=Y,-.^ix) 

lil<" 

For example sin(/i) = h — ^ ii h E so that h'^ = 0. 

It is possible to generalize several results of the present work to functions of class C" 
only, instead of smooth ones. However it is an explicit purpose of this work to simplify 
statements of results, definitions and notations, even if, as a result of this searching for 
simplicity, its applicability will only hold for a more restricted class of functions. Some more 



general results, stated for functions, but less simple can be found in Giordano 
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Note that m = f\x) G M, i.e. the slope is a standard real number, and that we can use 
the previous formula with standard real numbers x only, and not with a generic x G *M, but 



m 



we shall remove this limitation in subsequent works (see also Giordano 

If we apply this theorem to the smooth function p{r) := j^^"^ f(t) dt, for / smooth, then 
we immediately obtain the following result frequently used in several informal calculations: 

Corollary 23. Let A be open inM., x E A and f : A — > M smooth. Then 

fx+h 

yheD: / f{t)dt = h - fix). 

J X 

Moreover f{x) G M is uniquely determined by this equality. 



X. NILPOTENT INFINITESIMALS AND ORDER PROPERTIES 

Like in other disciplines, also in mathematics the layout of a work reflects the personal 
philosophical ideas of the authors. In particular the present work is based on the idea that 
a good mathematical theory is able to construct a good dialectic between formal properties, 
proved in the theory, and their informal interpretations. The dialectic has to be, as far as 
possible, in both directions: theorems proved in the theory should have a clear and useful 
intuitive interpretation and, on the other hand, the intuition corresponding to the theory 
has to be able to suggest true sentences, i.e. conjectures or sketch of proofs that can then 
be converted into rigorous proofs. 

In a theory of new numbers, like the present one about Fermat reals, the introduction 
of an order relation can be a hard test of the excellence of this dialectic between formal 
properties and their informal interpretations. Indeed if we introduce a new ring of numbers 
(like *M) extending the real field M, we want that the new order relation, defined on the 
new ring, will extend the standard one on M. This extension naturally leads to the wish of 
findings a geometrical representation of the new numbers, according to the above principle 
of having a good formal/informal dialectic. 

We want to start this section showing that in our setting there is a strong connection 
between some order properties and some algebraic properties. In particular, we will show 
that it is not possible to have good order properties and at the same time a uniqueness 
without limitations in the derivation formula. In the following theorem we can see that the 
property h ■ k = is a general consequence if we suppose to have a total order on D. 
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Theorem 24. Let {R, <) be a generic ordered ring and D C R a subset of this ring, such 
that 

1. Oe D 

2. yheD : h"^ = and -he D 

3. {D, <) is a total order 
then h ■ k = for every h, k e D. 

This theorem implies that if we want a total order in our theory of infinitesimal numbers, and 
if in this theory we consider D = {hlh"^ = 0}, then we must accept that the product of any 
two elements of D must be zero. For example, if we think that a geometric representation 
of infinitesimals is not possible if we do not have, at least, the trichotomy law, then in this 
theory we must also have that the product of two first order infinitesimals is zero. 
Proof: Let h, k E D he two elements of the subset D. By hypotheses 0, —h, —k G D, 
hence all these elements are comparable with respect to the order relation <, because, by 
hypotheses this relation is total in D. E.g. 

h < k or k < h 

We will consider only the case h < k, because analogously we can deal with the case k < h, 
simply exchanging everywhere h with k and vice versa. 

First sub-case: k > 0. By multiplying both sides of /i < A; by A; > we obtain 

hk < k'^ (23) 

li h > then, multiplying by > we have < hk, so from (l23l) we have < hk < k"^ = 0, 
and hence hk = 0. 

If /i < then, multiplying by /c > we have 

hk<0 (24) 

If, furthermore, h > —k, then multiplying by > we have hk > —k"^, hence form ( l24l l 
0>hk>-k^ = 0, hence hk = 0. 

If, otherwise, h < —k, then multiplying by — /i > we have —h'^ = 0<hk<0 from (l24l) . 
hence hk = 0. This concludes the discussion of the case k > 0. 
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Second sub-case: k < 0. In this case we have h < k < 0. Multiplying both inequalities by 
h <0 we obtain h'^ = > hk > and hence hk = 0. 

■ 

So, the trichotomy law is incompatible with the uniqueness in a possible derivation for- 
mula like 

3\meR: yheD: f{h) = /(O) + h-m (25) 

framed in the ring R of Theorem [241 In fact, if a, b E D are two elements of the subset 
D C then both a and b play the role of m G -R in (l25ll for the linear function 

f : h e D ^ h- a = e R 

So, if the derivation formula (l25ll applies to linear functions (or less, to constant functions), 
the uniqueness part of this formula cannot hold in the ring R. 

In the next section we will introduce a natural and meaningful total order relation on 'M. 
Therefore, the previous Theorem [24] strongly motivate that for the ring of Fermat reals 'R 
we must have that the product of two first order infinitesimals must be zero and hence, that 
for the derivation formula in 'R the uniqueness cannot hold in its strongest form. Since we 
will also see that the order relation permits to have a geometric representation of Fermat 
reals, we can summarize the conclusions of this section saying that the uniqueness in the 
derivation formula is incompatible with a natural geometric interpretation of Fermat reals 
and hence with a good dialectic between formal properties and informal interpretations in 
this theory. 

XI. ORDER RELATION 

From the previous sections one can draw the conclusion that the ring of Fermat reals 
•R is essentially "the little-oh" calculus. But, on the other hand the Fermat reals give us 
more flexibility than this calculus: working with 'R we do not have to bother ourselves 
with remainders made of "little-oh", but we can neglect them and use the powerful algebraic 
calculus with nilpotent inflnitesimals. But thinking the elements of 'R as new numbers, and 
not simply as "little-oh functions", permits to treat them in a different and new way, for 
example to deflne on them an order relation with a clear geometrical interpretation. 
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First of all, let us introduce the useful notation 

V°t > : V{t) 

and we will read the quantifier V°t > saying "for every t > (sufficiently) small", to 
indicate that the property V{t) is true for all t in some right neighborhood of t = (recall 
that, by Definition [2l our little-oh polynomials are always defined on ]R>o), i.e. 

3(5 > : Vt G [0,(5) : V{t). 

The first heuristic idea to define an order relation is the following 

X <y <^==^ X — y <Q <^=^ : 2; = in *M and x — y < z 

More formally: 

Definition 25. Let x, y E 'M, then we say 

X < y 

iff we can find z G 'M such that 2; = in 'M and 

V°t > : xt<yt + zt 

Recall that z = in 'M is equivalent to Zt = o{t) for t ^ 0~^. It is immediate to see that 
we can equivalently define x < y if and only if we can find x' = x and y' = y in 'M. such 
that Xt < yt for every t sufficiently small. From this it also follows that the relation < is 
well defined on 'M, i.e. ii x' = x and y' = y in 'M and x < y, then x' < y' (recall that, 
to simplify the notations, we do not use equivalence classes as elements of 'M but directly 
little-oh polynomials). As usual we will use the notation x < y for x < y and x ^ y. 

Theorem 26. The relation < is an order, i.e. is reflexive, transitive and anti- symmetric; it 
extends the order relation ofM and with it (*]R, <) is an ordered ring. Finally the following 
sentences are equivalent: 

1. he Doo, i.e. h is an infinitesimal 

2. Vr G M>o : -r < h < r 
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Hence an infinitesimal can be thought of as a number with standard part zero, or as a 
number smaller than every standard positive real number and greater than every standard 
negative real number. 

Proof: We only prove the prove the property 

X <y and w >Q =^ x ■ w < y ■ w, 

the others being a simple consequence of our Definition [251 Let us suppose that 

xt<yt + zt V°t>0 (26) 
wt > z[ V°t > 

then Wt — z[ > for every t small and hence from fl26l) 

xt ■ {wt - z[) < yt ■ {wt - z[) + zt ■ {wt - z[) V°t > 

from which it follows 

Xt ■ Wt < yt ■ Wt + i-xtz[ - ytz[ + ZtWt - Ztz[) V^t > 

But —xz' — yz' + zw — zz' = in 'M because z = and z' = and hence the conclusion 
follows. 

■ 

Example. We have e.g. dt > and dt2 — 3 dt > because for t > sufficiently small 
t^^'^ > 3t and hence 

- 3t > V°t > 0. 

From examples like these ones we can guess that our little-oh polynomials are always locally 
comparable with respect to pointwise order relation, and this is the first step to prove that 
for our order relation the trichotomy law holds. In the following statement we will use the 
notation V^t > : 'P(t), that naturally means 

yH>0: t^O =^ V{t) 

where V{t) is a generic property depending on t. 

Lemma 27. Let x, y E 'M, then 
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1. °x < °y 



V°t > : x,< 



Xt < Vt 



2. If °x = °y, then 



(V°t > : xt<yt) or (V°t > : Xt>yt) or [x = y in 'M) 



Proof: 

(21) Let us suppose that °x < °y, then the continuous function t>0\-^yt — XtE'K assumes 
the value yo — Xq > hence is locally positive, i.e. 

V°t>0: xt<yt 

[B) Now let us suppose that °x = °y, and introduce a notation for the potential decompo- 
sitions of X and y (see Definition [7l). From the definition of equality in 'M, we can always 
write 

N 

xt = °x + ^ai- 1"' + zt Vt > 

i=l 
M 

yt = °y + J2^rt'' +wt Vt > o 

where x = °x + J2iLi c^i ' t'^' y = °y + Pj ■ are the potential decompositions of 
X and y (hence < < ctj+i < 1 and Q < (3j < < 1), whereas w and z are little-oh 
polynomials such that Zt = o(t) and Wt = o(t) for t ^ 0^. 

Case: ai < hi In this case the least power in the two decompositions is ai ■ and 
hence we expect that the second alternative of the conclusion is the true one if ai > 0, 
otherwise the first alternative will be the true one if ai < (recall that always 7^ in 
a decomposition). Indeed, let us analyze, for t > 0, the condition Xt < yt'- the following 
formulae are all equivalent to it 

N N 

Y^a^■t''^ f^j ■ t^' +Wt-Zt 



i=l 



0=1 



±ai 



N 



ai 



ai 



1=2 
N 



< t 



ai 



N 



0.1 



U=i 



N 



1=2 j=l 
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Therefore, let us consider the function 

N N 
j=l i=2 

We can write 

-ai 



[Wt - Zt) -t ' = t 



and 



as i — > 0+ because Wt — o{t) and Zt — o{t). Furthermore, ai < 1 hence t^~°-^ 
is bounded in a right neighborhood of t — 0. Therefore, {wt — Zt) ■t~°-^ and the function 
/ is continuous at t = too, because Oj < Oj and ai < 6i < bj. By continuity, the function 
/ is locally strictly positive if and only if /(O) = —ai > 0, hence 

(V°t > : xt<yt) ^ ai < 
(V°t > : xt>yt) ^ ai > 

Case: ai > bi We can argue in an analogous way with bi and /3i instead of oi and cti. 
Case: Oi = bi We shall exploit the same idea used above and analyze the condition xt < yt- 
The following are equivalent ways to express this condition 



t 



ai 



N 



=2 

AT 



< t' 



N 



3=2 



N 



i=2 j=2 

Hence, exactly as we have demonstrated above, we can state that 



ai < Pi 
ai > Pi 



V°t>0: xt<yt 
V°t>0: xt>yt 



Otherwise ai — Pi and we can restart with the same reasoning using a2, &2, Q;2, P2, etc. If 
N — M, the number of addends in the decompositions, using this procedure we can prove 
that 

yt>0: xt^yt + wt- Zt, 

that is a; = y in 'M. 

It remains to consider the case, e.g., < M. In this hypotheses, using the previous 
procedure we would arrive at the following analysis of the condition Xt < yt- 

< X • +wt-zt 

j>N 



30 



< t^^+i 



Hence 



j>N+l 
j>N+l 

> ^ yH>0: xt<yt 

Pn+1 < =^ V't>0: xt>yt 



This lemma can be used to find an equivalent formulation of the order relation. 
Theorem 28. Let x, y E 'M, then 

1. X <y <^==^ (V°t > : xt < yt) or (x = y in *R) 

2. X <y <^=^ (V°t > : xt < yt) and {x ^ y m 'R) 
Proof: 

d) =^ If < °y then, from the previous Lemma [27l we can derive that the first alternative 
is true. If °x = °y, then from Lemma [27] we have 

(V°t >0: xt<yt) or {x = y in 'M) or (V°t > : Xt > yt) (27) 

In the first two cases we have the conclusion. In the third case, from x < y we obtain 

V°t > : xt<yt + zt (28) 

with Zt = o(t). Hence from the third alternative of ( l27l ) we have 

0<xt-yt<zt V°t > 

and hence lim(^o+ ~ 3; = y in 'M. 

[21 ) <^= This follows immediately from the refiexive property of < or from the Definition [25l 
[B) =^ From X < y we have x < y and x ^ y, so the conclusion follows from the previous 

m 

IE) From > : Xt < yt and from [11 it follows x < y and hence x < y from the 
hypotheses x ^ y. 

m 

Now we can prove that our order is total 
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Corollary 29. Let x, y E 'R, then in 'M we have 

1. X <y or y < X or x = y 

2. X < y or y < x or x = y 
Proof: 

(21) If °x < °y, then from Lemma [27] we have Xt < yt for t > sufficiently small. Hence 
from Theorem [28] we have x < y. We can argue in the same way if °x > °y. Also the case 
°x = °y can be handled in the same way using [B of Lemma [271 
[B) This part is a general consequence of the previous one. 

■ 

From the proof of Lemma [27] and from Theorem [28] we can deduce the following 
Theorem 30. Let x, y e *R. If°x ^ °y, then 

X < y <^==^ X < y 

Otherwise, if °x = °y, then 

1. If uj{x) > uj{y), then x > y ijf°Xi > 

2. Ifu{x) = uj{y), then 

°xi > °yi =^ x> y 
°xi < °yi =^ X <y 

Example. The previous Theorem gives an effective criterion to decide whether x < y oi 
not. Indeed, if the two standard parts are different, then the order relation can be decided 
on the basis of these standard parts only. E.g. 2 + dt2 > 3 dt and 1 + dt2 < 3 + dt. 

Otherwise, if the standard parts are equal, we firstly have to look at the order and at the 
first standard parts, i.e. °xi and °yi, which are the coefficients of the biggest infinitesimals 
in the decompositions of x and y. E.g. 3dt2 > 5dt, and dt2 > adt for every a G M, and 
dt < dt2 < dta < . . . < dtk for every k > 3, and dt^ > 0. 

If the orders are equal we have to compare the first standard parts. E.g. 3 dts > 2 dt^. 
The other cases fall within the previous ones, because of the properties of the ordered ring 
•M. E.g. we have that dtg - 2dt3 + 3dt < dig - 2dt3 + dtg/z if and only if 3dt < dtg/a, 
which is true because uj{ dt) = 1 < uj{ dt3/2) = |. Finally dts — 2 dta + 3 dt > dts — 2 dts — dt 
because 3 dt > — dt. 
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XII. ABSOLUTE VALUE, POWERS AND LOGARITHMS 

Having a total order we can define the absolute value in the usual way, and, exactly like 
for the real field M, we can prove the usual properties of the absolute value. Moreover, also 
the following cancellation law is provable. 

Theorem 31. Let he'M\ {0} and r, s eR, then 

\h\ ■ r < \h\ ■ s =^ r < s 

Proof: In fact if \h\ ■ r < \h\ ■ s then from Theorem [28] we obtain that either 

V°t > : \ht\ -r < \ht\ -s (29) 

or \h\ ■ r = \h\ ■ s. But h ^ so 

(V°t >0: ht>0) or (V°t > : ht < O) 

hence we can always find a t > such that \ht\ ^ and to which (1291 ) is applicable. 
Therefore, in the first case we must have r < s. In the second one we have 

\h\ ■ r = \h\ ■ s 

but h ^ 0, hence |/;,| 7^ and so the conclusion follows from Theorem [TTl 

■ 

Due to the presence of nilpotent elements in 'M, we cannot define powers and loga- 
rithms log^. y without any limitation. E.g. we cannot define the square root having the usual 
properties, like 

X e'R =^ Vx e'R (30) 
x = y m*R =^ = 7y in "M (31) 

Vx^ = \x\ 

because they are incompatible with the existence oi h E D such that /i^ = 0, but h 0. 
Indeed, the general property stated in the Subsection [TV] permits to obtain a property like 
( l30l) (i.e. the closure of *R with respect to a given operation) only for smooth functions. 
Moreover, the Definition [I9] states that to obtain a well defined operation we need a locally 
Lipschitz function. For these reasons, we will limit to x > and x invertible only, and 
log^y to X, y > and both x, y invertible. 
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Definition 32. Let x, y E *M, with x strictly positive and invertible, then 

1. xy := [t>0^xfUi^.^ 

2. If y > and y is invertible, then log^ y ■= [t > ^ logxtyt\= in •» 
Because of Theorem [28] from x > we have 

V°t > : > 

so that, exactly as we proved in Subsection [IV] and in Definition [T9l the previous operations 
are well defined in 'M because °x 7^ 7^ °y. From the elementary transfer theorem [20] 
the usual properties follow. To prove the usual monotonicity properties, it suffices to use 
Theorem [M 

Finally, it can be useful to state here the elementary transfer theorem for inequalities, 
whose proof follows immediately from the definition of < and from Theorem [28l 

Theorem 33. Let A he an open subset of M", and t, a : A — > M be smooth functions. 
Then 

yxe'A: 't{x) < 'a{x) 

iff 

Vr G A : r(r) < cr(r). 

XIII. GEOMETRICAL REPRESENTATION OF FERMAT REALS 

At the beginning of this article we argued that one of the conducting idea in the con- 
struction of Fermat reals is to maintain always a clear intuitive meaning. More precisely, 
we always tried, and we will always try, to keep a good dialectic between provable formal 
properties and their intuitive meaning. In this direction we can see the possibility to find a 
geometrical representation of Fermat reals. 

The idea is that to any Fermat real x G "M we can associate the function 

N 

teR>o^°x + J2 ■ t^^'^'^''^ e ^ (32) 
1=1 

where is, of course, the number of addends in the decomposition of x. Therefore, a 
geometric representation of this function is also a geometric representation of the number 
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x, because different Fermat reals fiave different decompositions, see Tfieorem[5l Finally, we 
can guess that, because the notion of equality in 'M depends only on the germ generated 
by each little-oh polynomial (see Definition [4j), we can represent each a; G "M with only the 
first small part of the function (i32l ). 

Definition 34. If x G 'M and 6 G M>o, then 

f 

graphsix) := U°x + Yl ' t^^'^'^^^t) \0<t<6 
[ i=i 

where is the number of addends in the decomposition of x. 

Note that the value of the function are placed in the abscissa position, so that the correct 
representation of graph^(x) is given by the figured 



t 50 - 




2 4 6 8 10 12 14 16 18 20 

FIG. 1: The function representing the Fermat real dt2 G D-^ 

This inversion of abscissa and ordinate in the graph^(a;) permits to represent this graph as 
a line tangent to the classical straight line R and hence to have a better graphical picture. 
Finally, note that if x G M is a standard real, then N = and the graph^(x) is a vertical 
line passing through °x = x. 

The following theorem permits to represent geometrically the Fermat reals 
Theorem 35. If S E M>o, then the function 

X graphs (x) C 
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is injective. Moreover if x, y G 'M, then we can find 5 G M>o ( depending on x and y ) such 
that 

X < y 

if and only if 

yp,q,t: G graphs (x) , e graphs (y) =^ p<q (33) 

Proof: The application p(x) := graph5(x) for x G "M is well defined because it depends 
on the terms °x, °Xi and uji{x) of the decomposition of x (see Theorem [5] and Definition [8]). 
Now, suppose that graph5(a;) = graphed/), then 

N M 

Vt G [0, 6): °x + Y, ■ i^'^'^''^ = °y + Y^ °yj ■ t^'"''^''^. (34) 

i=l j=l 
Let us consider the Fermat reals generated by these functions, i.e. 



X 



N 



t>Q^°X + ^°Xi-t 



i=l 
AI 



y:= t>o^°y + J2°yrt' 

then the decompositions of x' and y' are exactly the decompositions of x and y 



N 



x' = °x 



X 



i=l 
M 



(35) 
(36) 



But from (l34ll it follows x' = y' in 'M, and hence also x = y from (l35ll and (l36ll . 

Now suppose that x < y, then, using the same notations of the previous part of this 
proof, we have also x' = x and y' = y and hence 

N M 

o I 
X = X + 

i=l j=l 

We apply Theorem [28l obtaining that locally x[ < y[, i.e. 

N M 

35 > : V°t > : °x + ^°Xi- 1^/^'^^^ <°y + Y^ °yj " t^'^'^^\ 

i=i j=i 

This is an equivalent formulation of (l33l l. and, because of Theorem [281 it is equivalent to 

x' = X < y' = y. 
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Example. In figure [2] we have the representation of some first order infinitesimals. 




FIG. 2: Some first order infinitesimals 

The arrows are justified by the fact that the representing function fl32l) is defined on R>o 
and hence has a clear first point and a direction. The smaller is a G (0, 1) and the nearer 
is the representation of the product a dt, to the vertical line passing through zero, which is 
the representation of the standard real x = 0. Finally, recall that dtfc ^ D \i and only if 
1 < < 2. 

If we multiply two infinitesimals we obtain a smaller number, hence one whose representation 
is nearer to the vertical line passing through zero, as represented in figure [3] 



dt4 • dt3 




FIG. 3: The product of two infinitesimals 

In figured] we have a representation of some infinitesimals of order greater than 1. We can 
see that the greater is the infinitesimal h E Da (with respect to the order relation < defined 
in 'R) and the higher is the order of intersection of the corresponding line graph^(/i). 

Finally, in figure [5] we represent the order relation on the basis of Theorem [35l 
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FIG. 5: Different cases in which Xj < yi 

Intuitively, the method to see if x < y is to look at a suitably small neighborhood (i.e. at a 
suitably small 5 > 0) at t = of their representing lines gTaph^{x) and gTaph^{y): if, with 
respect to the horizontal directed straight line, the curve graph^^x) comes before the curve 
graphed/), then x is less than y. 



XIV. SOME ELEMENTARY EXAMPLES 



The elementary examples presented in this section want to show, in a few rows, the 
simplicity of the algebraic calculus of nilpotent infinitesimals. Here "simplicity" means that 
the dialectic with the corresponding informal calculations, used e.g. in engineering or in 
physics, is really faithful. The importance of this dialectic can be glimpsed both as a proof 
of the fiexibility of the new language, but also for researches in artificial intelligence like 
automatic differentiation theories (see e.g. Griewank and references therein). Last but 
not least, it may also be important for didactic or historical researches. Several examples are 
directly taken from analogous of Bell and the reader is strongly invited to compare the 
two theories in these cases. In particular, in our point of view, is not positive, like in some 
parts of Bell Q], to return back to a non rigorous use of infinitesimals. Mathematical theories 
of infinitesimals, like our ring of Fermat reals or NSA or SIA, are great opportunities to avoid 
several fallacies of the informal approach (our discussion in Section |X] is a clear example), and 
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to advance further, with the new knowledge originating from the rigorous theory, opening 
the possibihty to use infinitesimal methods in more general, and less intuitive, frameworks 
(like e.g. infinite dimensional spaces of mappings, see Giordano jl^). Once again, the key 
point is the dialectic between formal and informal thoughts and not a single part only. 



A. The heat equation. 

In this and the following section we simply use the language of 'M to reformulate the cor- 



responding deductions of Vladimirov [28|. Let us consider a body (identified with its local- 
ization) B C M.^ and denote with Ib '■= int(i?) its interior. On Ib are given the smooth 
functions p : Ib — > M, c : — > K and k : Ib — > M, interpreted respectively as the mass 
density, the specific heat capacity and the coefficient of thermal conductivity. Let us note 
that assuming these functions as defined on Ib without any favored direction corresponds 
physically to assume that B is an isotropous body. Moreover, let u : Ib ^ [0, +00) — > M be 
the smooth function representing the temperature of the body B at each point x E Ib and 
time t G [0, +00). To deduce the heat diffusion equation we fix an internal point x E Ib and 
an infinitesimal volume V . More precisely, we say that a subset of 'M^ of the form 

V = V{x,5x) = {y -5xi<2{y-x)-ei<5xi Vi = l,2,3} (37) 

is an infinitesimal parallelepiped liSv := 5xi-5x2-Sx^ G -Dqo, i-e. if the corresponding volume 
is an infinitesimal of some order. Here (ci, 62, 63) is the natural base of and notations of 
the form 5y G *M are only useful to underline that the infinitesimal increment is associated 
to the variable y: here 5 is not an operator and we use it instead of the common dy to avoid 
confusion with our dy introduced in Definition [3 Because x E Ib, the inclusion V C'B 
follows, so that V can be thought as the sub-body of B corresponding to the infinitesimal 
parallelepiped parallel to coordinate axis and centered at x. This sub-body V interacts 
thermally with its complement CV := 'B \V and with external sources of heat. In the 
infinitesimal time interval 6t G -Doo, the sub-body V exchanges with its complement CV the 
heat flowing perpendicularly to the surface of V (Fourier's law): 
3 



Qcvy = St ■ ^ 



i=l 



k(x + Shi) ■ 7— (2^ + Shi, t) — k(x — Shi) ■ 7— (2^ — Shi, t) 

OCi OCi 



(38) 
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where Shi := ^Sxt ■ Cj G 'M and 5si := lljyi^^i ^ Choosing the infinitesimals so that 

6v-6te D, 

we have that 6t ■ 6si ■ {6xi)'^ = 6t ■ 6v ■ Sxi = from Theorem [12] (e.g. we can choose 
6xi = dte and 6t = dt2). From this and the use of infinitesimal Taylor's formula in (l38l l. 
simple calculations give 

Qcv,v = div [k ■ grad(M)] (x, t) ■ 6v ■ 6t. (39) 

Of course, this calculations correspond to the infinitesimal version of the Gauss-Ostrogradskij 
theorem. Interacting thermally with external sources, the sub-body V exchanges the heat 

Qe.t,v = F{x,t) ■ 5v 5t, (40) 

where F : Ib y< [0, +oo) — > R is a smooth function representing the intensity of the thermal 
sources. The total heat Qcv,v + Qexty corresponds to an increasing of temperature of V 
equal to u{x, t + 6t) — u{x, t) and hence to an exchange of heat with the environment equal 
to 

Qenv,V = [u{x, t + 6t) - u{x, t)] ■ c(x) ■ p(x) ■ 6v = QcV,V + Qext,V- (41) 

From this and ( l39ll . (I40l ). the infinitesimal Taylor's formula and the cancellation law we 
obtain the conclusion: 

du 

c{x) ■ p{x) ■ t) = div [k ■ grad(-u)] (x, t) + F{x, t). 

To stress that the previous deduction is now completely rigorous we can now state the 
following theorem, without any mention to the physical interpretation: 

Theorem 36. Let B C M.'^ and Ib '■= int(i?) its interior. Let us consider the smooth 
functions p : Ib — ^ R, c : Ib — > M., k : Ib — ^ m : x [0, +00) — > M and 
F : Ib ^ [0, +00) — > M. Finally let us consider a point {x,t) G Ib ^ [0, +00) and define V, 
Qcvy , Qexty , Qerwy (IS in (l37l ). (l38l) . (l40l) and (1411) . where 6v ■ 6t E D. Then it results 

Qcvy + Qext,y 

if and only if the following relation holds 

du 

c{x) ■ p{x) ■ -^{x, t) = div [k ■ grad(M)] {x, t) + F{x, t). 
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Unfortunately, this statement does not sufficiently underline the great difference that 
takes place between the physical content in the definition of Qcvy^ i-e. the Fourier's law, 
and that in the definition of Qexty- In an axiomatic framework for thermodynamics (see e.g. 
Truesdell j^l), the notion of heat flux Qab going from a body A to a body B can be taken 
as primitive; in that case ( !38ll becomes an important assumption, whereas ( l40l) is simply the 
deflnition of the intensity F{x,t) = 



In elementary Physics, an electric dipole is usually defined as "a pair of charges with op- 
posite sign placed at a distance d very less than the distance r from the observer". Conditions 
like r ^ d are frequently used in Physics and very often we obtain a correct formalization 
if we ask G *M infinitesimal but r G M \ {0}, i.e. r finite. Thus we can define an electric 
dipole as a pair (^1,^2) of electric particles, with charges of equal intensity but with opposite 
sign such that their mutual distance at every time t is a first order infinitesimal: 



In this way we can calculate the potential at the point x using the properties of D and using 
the hypothesis that r is finite and not zero. In fact we have 



B. Electric dipole. 



Vt: \p,{t) - p2it)\ =: \dt\ =:dteD. 



(42) 




Ti : = 



and if f := r2 — f then 




because for (l42l) d'^ 



For our hypotheses on d and r we have that 



r • d 



G D hence from 







the derivation formula 




In the same way we can proceed for 1/ri, hence: 




q r • d 



4:7160 
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The property (P — is also used in the calculus of the electric field and for the moment of 
momentum. 

C. Newtonian limit in Relativity. 

Another example in which we can formalize a condition like r ^ d using the previous 
ideas is the Newtonian limit in Relativity; in it we can suppose to have 

• yt : VteD2 and c e R 

• Vx e M4 : gij{x) — Tjij + hij{x) with hij{x) e D. 

where is the matrix of the Minkowski's metric. This conditions can be interpreted as 

Vt ^ c and hij{x) <^ 1 (low speed with respect to the speed of light and weak gravitational 
field). In this way we have, e.g. the equalities: 

—== = 1 + — and ^/l^^h^{x) ^1- -h4i{x). 

D. Linear differential equations. 

Let 

d^T/ dy 
L{y) := + . . . + ^^-1^ +A^-y^Q 

be a linear differential equation with constant coefficients. Once again we want to discover 
independent solutions in case the characteristic polynomial has multiple roots e.g. 

(r — r^Y ■ (r — Tg) ■ . . . ■ (r — rjv) = 0. 

The idea is that in 'R we have {r-rif = also if r = ri+/i with heD. Thus y{t) = e(^i+'*)* 
is a solution too. But e^''^+''^* = e*"!* + ht ■ e^^*, hence 

L = 

= L [e^i* + ht ■ e''^*] 

= L [e'''] +h-L[t- e'^'] 

We obtain L • e^^*] = 0, that is yi{t) — t ■ e^^* must be a solution. Using k-ih order 
infinitesimals we can deal with other multiple roots in a similar way. 
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E. Circle of curvature. 

A simple application of the infinitesimal Taylor's formula is the parametric equation 
for the circle of curvature, that is the circle with second order osculation with a curve 
7 : [0, 1] — > M'^. In fact if r G (0, 1) and 7^ is a unit vector, from the second order 
infinitesimal Taylor's formula we have 

h'^ h? ^ 

Wh e D2 : 7(r + h) = 7^ + /l7r + — 7r = 7r + ^^r + —Cr fir (43) 

where n is the unit normal vector, t is the tangent one and Cr the curvature. But once 

2 1,2 

again from Taylor's formula we have sin(c/i) = ch and cos(c/i) = 1 — Now it suffices to 
substitute h and ^ from these formulas into fl43l) to obtain the conclusion 

\/h ^ D2 : 7(r + h) = f 7^ H — -] H ■ \?,m.icrh)tj. — co?,iCrh)nr~\ . 

V CrJ Cr 

In a similar way we can prove that any / G C°°(M, M) can be written Wh G as 

k k 

f{h) = a„ ■ cos{nh) + 6„ ■ sm{nh), 

n=0 n=0 

so that now the idea of the Fourier series comes out in a natural way. 

F. Commutation of differentiation and integration. 

n n 

This example derives from Kock [2C!||, Lavendhomme [22]. Suppose we want to discover the 
derivative of the function 

9[x 



:= / f{x,t)dt Vx G 

J a(x) 



where a, [3 and / are smooth functions. We can see (7 as a composition of smooth functions, 
hence we can apply the derivation formula, i.e. Theorem | 



g{x + h) = / f{x + h,t)dt 

Ja(x+h] 



a{x+h) 

a{x) na(x) Qr 

f{x,t)dt + h- / -^{x,t)dt+ 

a{x)+ha'{x) J a(x)+ha'{x) C)X 



nl3{x) nt^ix) 

+ / f{x,t)dt + h- / 

J a(x) J a 



(x, t) dt+ 



J ci(x) c)x 

P{x)+hf3'{x) i.p{x)+hf3'{x) Qn 

fix,t)dt + h- / -^1 

l3{x) Jf3(x) OX 
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Now we use /i^ = to obtain e.g. (see Corollary [23l) : 



h- / —(x,t)dt = —h'^-a'{x)-—{a{x),t) = 

Ja{x)+ha'{x) 



f 

Ja( 



and 

f{x, t) dt = —h ■ a'{x) ■ f{a{x),t). 

' a{x)+ha' (x) 

Calculating in an analogous way similar terms we finally obtain the well known conclusion. 
Note that the final formula comes out by itself so that we have "discovered" it and not simply 
we have proved it. From the point of view of artificial intelligence or from the didactic point 
of view, surely this discovering is not a trivial result. 

G. Schwarz's theorem. 

Using nilpotent infinitesimals we can obtain a simple and meaningful proof of Schwarz's 
theorem. This simple example aims to show how to manage some differences between our 
setting and SDG. Let / : V — ^ E he aC^ function between spaces of type V = M™, E = W^ 
and a E V, we want to prove that d^/(a) : V x V — > E is symmetric. Take 

h, j e Poo 

jkh E D-LQ 

(e.g. we can take kt = dt2,ht = jt = dt4 so that jkh = dt, see also Theorem [T2ll. Using 
k E D2, we have 



j ■ f{x + hu + kv) 



3 ■ 



k^ 

f{x + hu) + k dyf{x + hu) + —dlf{x + hu) 



(44) 



= 3 ■ + ^"w) + 3^ ■ dvf{x + hu) 
where we used the fact that k"^ E D and j infinitesimal imply jk^ = 0. Now we consider 
that jkh E D so that any product of type jkhi is zero for every i E Doc, so we obtain 

jk ■ dyf{x + hu) = jk ■ d^f{x) + jkh ■ du{d^f){x). (45) 

But k E D2 and jk"^ = hence 

j ■ f{x + kv) - j ■ f{x) = jk ■ d^f{x). 
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Substituting this in (14511 and hence in (1441) we obtain 

j ■ [f{^ + hu + kv) — f{x + hu) — f{x + kv) + /(x)] = 

(46) 

= jkh ■ du{d^f){x). 

The left hand side of this equahty is symmetric in u, v, hence changing them we have 

jkh ■ du{d^f){x) = jkh ■ (9„(9„/)(x) 

and thus we obtain the conclusion because jkh ^ and du{dvf){x), dv{duf){x) G E. From 
(l46ll it follows directly the classical limit relation 

f{x + htu + ktv)-f{x + htu)-f{x + ktv) + f{x) r. r. . 

hm — = dudyf{x) 

t^o+ htkt 

H. Area of the circle and volumes of revolution. 

A more or less meaningful proof of the familiar formula for the area of a circle depends 
on what axioms are assumed and how much general the definitions are. In this example 
we want to show the possibility to define suitable smooth functions using an infinitesimal 
property. Let us assume the axioms for the real field R; prove from them the existence of 
the smooth functions sin and cos; define vr as a suitable zero of these functions (see e.g. 
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26|) and define the length of an arc of circle of radius r, parametrized by 



x{6) = r ■ cos{9) and y{6) = r ■ sin(6'), as the unique function s that verifies 

[s{e + k) - s{e)f = [x{e + k) - x{e)f + [y{e + k) - y{e)f v^gmv^gDs (47) 

s(0) = 0. (48) 

This definition can be justified in the usual way using a (second order!) infinitesimal right- 
angled triangle. The uniqueness of s follows from ( l47l ) and (l48l) . the smoothness of a; and y, 
the second order infinitesimal Taylor's formula and the cancellation law (Theorem [TTll : 

e ■ s{e) = x{e) ■ e + y{e) \/ke D2. 

From this and (i48ll we obtain the usual formula for s that, in our particular case, gives 
s{6) = r-6. Now we can think the area A{6 + h) — A{6) of a first order infinitesimal sector of 
the circle as the area of the isosceles triangle with sides of length r and base s{9 + h) — s{9). 
In fact, if P{9) = {rsm9,rcos9), then P{9 + h) = P{9) + h ■ t{9), where Tis the tangent 
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vector, so that in [9,9 + h], h ^ D, the circle is made of hnear segments. Therefore, the area 
74(6*) can be defined as the unique function that verifies 

A{9 + h)- A{9) = ^ [s{9 + h)- s{9)] ■ r cos \/9ER\/heD 
A{0) = 0. 

From this and the derivation formula we get 

h ■ A'{9) = hir ■ s\9) 
1 

A{9) = 2 / ' ^i^) 







In our case we get A{9) = ^r^ ■ 9 and hence the searched formula for 9 = 2tt. 
Analogously we can prove the familiar formula for volumes of revolution of a parametrized 
curve 7(m) = {x{u),y{u)), u G [a,h], around the x-axis. Let us define the volume as the 
unique smooth function V that verifies 

V{u + h)- V{u) = h-n- y{uf + ]^[h ■ ^ ■ y{u + hf - h ■ ^ ■ y{uf] (49) 

y(0) = (50) 

for every u G [a, 6] and h E D. This definition can be intuitively justified saying that the 
volume of the sector of revolution between u and u + h can be calculated as the sum of the 
cylinder of radius y{u) and height h plus one half of the difference between the cylinder of 
radius y{u + h) and height h and that of radius radius y{u) and the same height. Implicitly, 
we are using the straightness of the curve 7 in [u, u + h]. From (l49l) and the property h? = Q 
we easily obtain that V'{u) = vr • y{uy and hence the usual formula using (l50ll . 

I. Curvature. 

Let us consider the usual smooth parametrized curve 7(m) = {x{u),y{u)) for u G [a,b]. 
Let ip{u) G [0, vr] be the non-oriented angle (i.e. the one defined by the scalar product) 
between the tangent vector t = (i, y) and the unit vector i of the x-axis, so that 

\/ x^ + y"^ ■ cos if = X. 

Multiplying this equality by sin (p we easily obtain 

y ■ cos ip = X ■ simp. (51) 
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It is well known that the curvature of 7 at the point u G [a, b] can be calculated as the 
rate of change of the non-oriented angle (p{u) with respect to an infinitesimal variation in 
arc length s{u) defined by the analogous of (l47ll and f l48l) . These "rate of changes" can be 
defined in 'M as the unique (if it exists) standard c{u) G M defined by 

c{u) ■ [s{u + h) — s{u)] = (p{u + h) — (p{u) \/h G D. 

Indeed, from the cancellation law, i.e. Theorem \17\ there exists at most one such c{u) G M 
verifying this property. Because of this uniqueness we can also use the notation 

c{u) = + (52) 
^ ^ s{u + h) - s{u) ^ ' 



These ratios generalize the usual ratios between real numbers (see Giordano [l^ for more 
details). From f l52l ) and the derivation formula we get c{u) = -l^f?^ = whatever 
h G -D^o we choose. From this and the relation (l5TD (without using infinitesimals, but using 
standard differential calculus) we can obtain the usual formula c = ^^2^2^/2 ^it each point 
u G [a, b] where ip{u) 7^ | and •y^u) 7^ 0. 



J. Stretching of a spring (and center of pressure). 

If / : [a, b] — > M is a smooth function and we define J(x) := f{s) ds, then Corollary 
[23l and a trivial calculation with the derivation formula give 

J{x + h)- J{x) = i [fix + h) + fix)] V/i G D. (53) 

The right-hand side of (l53l ) is interpreted as the average value of / in the infinitesimal 
interval [x, x + h]. Analogous equalities can be obtain in the (i-dimensional case using 
suitable generalizations of the above cited corollary: e.g. if c? = 2 we have to use 

/ / fix,y)dxdy = hk-fi<d,Q) \/h, k E : h ■ k e D. 
Jo Jo 

These equalities are used by Bell |3] to calculate the center of pressure of a plane area and the 
work done in stretching a spring. The meaningfulness of such examples is however doubtful 
because they can be summarized saying: assume to have a smooth J satisfying ( l53l l: deduce 
from this and from the assumption J(0) = that J'(x) = fix). There is no real use of 
infinitesimals in this type of reasoning in every case where the definition J(x) := /(s) ds 
is customary, like in the cited examples. 
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K. The wave equation. 
The deduction of the wave equation in the framework of Fermat reals is very interesting 



for two main reasons. Firstly, in the classical deduction (see e.g. Vladimirov [28||) there 
are some approximations tied with Hook's law. Is it possible to make them rigorous using 
*M? Do we gain something using this increased rigour? E.g.: how can we formalize the 
approximated equalities used in the classical deduction? In what a sense is the wave equation 
an approximated equality valid for small oscillations only? 

Secondly, at the end of our deduction we will stress the physical principles as important 
mathematical assumptions of a suitable theorem. We are hence naturally taken to ask if 
these natural assumptions (some of which formulated using the infinitesimals of 'M) really 
have a model. In this way, we will see that no standard smooth function can satisfy these 
hypothesis, but we are forced to consider a non-standard one. E.g. f{x) = h ■ sin(a;) for 
X G "M and h G -Doo is an example of a non-standard smooth function; let us note that it 
is obtained by the standard smooth function g{y,x) := y ■ sin(x), x, y G M, by extension to 
and fixing one of its variables to a non-standard parameter h G Doo- 

f{x) = 'g{h,x) VxG'M. 

This will motivate strongly the further development of the theory of Fermat reals, in the 
direction of a more general theory including also these new smooth non-standard functions. 

Let us start considering a string making small transversal oscillations around its equi- 
librium position located on the interval [a, b] of the x axis, for a, 6 G M, a < 6. By 
hypotheses, string's position St C 'M^ is always represented by the graph of a given curve 7 : 
[a, b] X [0, +00) — > 'R^ (where [a,b] = {x e *R \ a < x < b} and [0, +00) = {x G 'M | < x}; 
in the following, we will always use these notations for intervals to identify the corresponding 
subsets of 'M, and not of M, and we will also use the notation jxt '■= 7(2;, t)): 

st = {-f^te*R^\a<x <b} VtG[0,+oo). 

Moreover, the curve 7 is supposed to be injective with respect to the parameter x G (a, 6): 

7xit lx2t Vt G [0, +00) Vxi, X2 G (a, b) : Xi 7^ Xa, 



so that the order relation on (a, 6) implies an order relation on the support St. For every 
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pair of points p = jxpt, Q = Ix^t £ St on the string at time t, we can define the sub-bodies: 

~P ■= {ixt \xp<x<h} 
P" := {ixt \a<x<Xp} 
pq ■= {lxt\xp <x<Xg} 

corresponding respectively to the parts of the string that follows the point p E St, that 
precedes the same point and that lies between the point p E St and the point q E St. 
It is usually implicitly clear that e.g. every sub-body of the form J)' exerts a force on 
each sub-body with which it is in contact, i.e. of the form pq or *p. Moreover, the force 
F{A, B) E "M^ that the sub-body A exerts on the sub-body B verifies the following equalities 



(see e.g. Truesdell [271): 

F{pq,*P) =F{^,'p) (54) 

F("g ,m) = F( g", V) (55) 

F{*p ,pq) = —F{pq,*p ) (action-reaction principle) (56) 

for every pair of points p, q E St and every time t E [0, +oo). Using this formalism, the 

tension at the point 'jxt £ St at time t E [0, +oo) can now be defined in the following way 

Tix,t):=F{^t,l^t)- (57) 



Now, let us consider the infinitesimal sub-body x,x + Sx := "fxt'tx+Sx^t ^ St located at time 
t between the points •jxt £ St and 'yx+Sx,t G St, where 5a; G -D is a generic first order 
infinitesimal. On this infinitesimal sub-body, mass forces of linear density G : [a, b] x 
[0, +oo) — y act, so that Newton's law can be written as 



P- ^x- TTiT = Fi^xt, x,x + 6x) + F{'yx+5x,t, x,x + 6x) + G ■ p ■ 6x, (58) 

where p : [a, b] x [0, +oo) — > 'M is the linear mass density and where, if not otherwise 
indicated, all the functions are calculated at (x, t) E (a, b) x [0, +oo). Of course, the contact 
forces appearing in Newton's law are due to the interaction of the infinitesimal sub-body 
with other sub-bodies in contact with its border 



d 



x,x + Sx 



{ixtj'yx+sxA ^ •] 
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Using action-reaction principle (l56ll and the equality (l55ll . with q = 'jx+5x,t and p = ■jxt so 



that pq = x,x + 6x, from (l58l) we have 



p- Sx- — = -F(x, x + 6x, 7^t) + F(7a;+5a;,t, 'yx+Sx,t) + G ■ p ■ Sx. 
Using (l54ll and the definition ( l57l l of tension we get 

= -T(x,t) + T(x + 5a;,t) + G-p-5x. (59) 

Up to this point of the deduction we have not used neither the hypotheses of small oscillations 
nor that of transversal oscillations. The second one can be easily introduced with the 
hypotheses 

G{x,t)-ei = Vx,t, (60) 

where (61,62) are the axis unit vectors. Using the notation f{x,t) for the non-oriented 
angle between the tangent unit vector t{x,t) at the point '^xt and the x axes (see (15T1)). the 
hypotheses of small oscillations can be formalized with the assumption 

ip{x,t) e D Vx,t. (61) 

This will permit to reproduce the classical deduction in the most faithful way (even if, as we 
will see later, a weaker assumption can be considered). Moreover, in the classical deduction 
of the wave equation, one considers only curves of the form •jxt = {x,u{x,t)). In this way 
from (EH) and the derivation formula we have 

--— ■ cos (p = smip 
ox 



dx 

' 972 ^ 2 



if e D 



so that (gf-) = and hence the total length of the string becomes: 



b 



1 + 



dx = b-a VtG[0,+oo). (62) 



By Hook's law, this justifies that the tension can be assumed to have a constant modulus 
T, not depending neither by the position x nor by the time t: 

T{x,t) =T ■t{x,t) Vx e (a,6) Vt G [0,+cx)). (63) 
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A tension T parallel to the tangent vector is the second part of the hypothesis about non 
transversal oscillations of the string. Let us note explicitly that the only standard continuous 
function verifying the equality L = b — a is the constant one, so the function u : [a, b] x 
[0, +00) — > *]R has to be understood as a non-standard one; later we will do further 
considerations about this important point. Projecting the equation ( 1591 ) on the y axis, we 
obtain 

p ■ 6x ■ = ■ ^{^^ t) ■ 6*2 + T ■ t(x + 6x,t) ■ 62 + G ■ 62 ■ p ■ Sx 

= —T sin {p{x, t) + T ■ sin ip{x + Sx, t) + G ■ p ■ 5x. 

But sin (y9 = = 1^ because G -D is a first order infinitesimal, hence 



cP'u 

p.6x. — = T. 



du , ^ , du 
-ix + Sx,t)--{x,t) 



+ G ■ p-5x 



■ 6x. (64) 



We cannot use the cancellation law with 5a; G -D to obtain the final result, because, as we 
mentioned above, the function u{x, t) G *M can assume non standard values, so it is time to 
clarify some points. As mentioned above, there does not exist a standard smooth function 
verifying all the assumptions or the physical principles we have used. Of course, everything 
depends by how we formalize the classical informal deduction used in elementary physics: 
e.g. we have chosen to use an equality sign in ( !62ll instead of an approximated equality; 
anyway we have to consider that if we use ~ to write (l62l) . then the problem becomes how to 
make more precise, physically, numerically or mathematically, this approximation; moreover, 
if we use an approximation sign in (l62ll , then we consistently must use the same sign both in 
( l63l) and therefore in the final wave equation. Nevertheless, smooth non standard functions 
can verify all the hypothesis and physical principles we have considered: e.g. the function 
u{x,t) := Mosin(x + u ■ t) is one of these if the maximum amplitude Uq E D and if p is 
constant, G = and T = ou'^p. 

Definition 37. If X C 'W and Y C 'Ry then we say that 

/ : X > Y is (non standard) smooth 

iff / maps X in Y and for every Xq G X we can write 

f{x) = 'g{p,x) yxe'VnX (65) 
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for some 

V open in W such that xq E *V 
p G *U, where U is open in W 
g G C^{U X V,W), 

where (-, -) : {[x]^,[y]^) e 'U x 'V i — > [{x,y)]^ G *{U x V) (see Definition H for the 
relation ~). 

In other words locally a smooth function / : X — > Y from X C to F C *W is 
constructed in the following way: 

1. start with an ordinary standard function g G C°°{U x V, My), with U open in MP and V 
open in M^. The space has to be thought as a space of parameters for the function 
9; 

2. consider its Fermat extension obtaining *g : *{U x V) — > 'My-, 

3. consider the composition *g o (— , — ) : *f/ x *V — > 'M^, where (— , — ) is the isomor- 
phism •[/ X ~ '(U X V) defined by {[x]^, [y]^) = [{x,y)]^; we will always use the 
identification 'U x 'V = *{U x V), so we will write simply x) instead of x) . 

4. fix a parameter p G 'f/ as a first variable of the previous composition, i.e. consider 
'g{P-, ~) '■ — ^ Locally, the map / is of this form: / = 'g{p, — ) = 'g{p, —)■ 

Because p = °p + h, with h G -Dqo, applying the infinitesimal Taylor's formula to variable p 
for the function *g{p,x) it is not hard to prove the following Theorem, that clarifies further 
the form of these non standard smooth functions, because it states that they can be seen 
locally as "infinitesimal polynomials with smooth coefficients": 

Theorem 38. Let X C'W and f : X — ^ 'R'" a map. Then it results that 

f : X — > *R" is non standard smooth 
if and only if for every xq E X vje can write 

/(x) = ^ a<,(x) -p^ \/xe'VnX, (66) 

\q\<k 

for suitable: 



52 



1. d, keN 



2. peDf 



3. V open subset ofW such that Xq G 'V 
I {a,)\,\<k family ofC^{V,W). 

In other words, every smooth function / : X — > *W can be constructed locally starting 
from some "infinitesimal parameters" 

Pi,...,PdeDk 

and from ordinary smooth functions 

and using polynomial operation only with pi, pd and with coefficients aq{—). Roughly 
speaking, we can say that they are "infinitesimal polynomials with smooth coefficients. The 
polynomials variables act as parameters only". 

As it is natural to expect, several notions of differential and integral calculus, including 
their infinitesimal versions, can be extended to this type of new smooth function (for more 



details, see the preprint Giordano ll6|), and these results will be presented in future works. 
In this sense, this deduction of the wave equation motivates strongly the future development 
of the theory of Fermat reals. 

On the other hand, we have to understand what type of cancellation law we can apply to 
fl64l) . For this end, we have to define the notion of equality up to fc-th order infinitesimals: 

Definition 39. Let m = °m + Yl!i=i ' ^'tuj,{m) be the decomposition of m G 'M and 
k G R>o U {oo}, then 

N 



Lkm := Lk{m) := °m + ^ °mi ■ dt<^^(m). 



Li}i(m)>k 



Finally if x, ?/ G *M, we will say x =fc |/ iff LkX = LkU in 'M, and we will read it as x is 
equal to y up to k-th order infinitesimals. 
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In other words, as it is easy to prove, we have 

X =k y °x = °y and u!{x — y) < k. 

Therefore, if we denote with 

h ■■= {x e Doo I uj{x) < k} , 

the set of all the infinitesimal of order less that or equal k (let us note that Ik C Dk), then 
we have that x =k y if and only if x — y G /fc. Equality up to k-th order infinitesimal is of 
course an equivalence relation and preserves all the ring operations of 'M. More in general 
these equalities are preserved by smooth functions / : *M — > *]R: 

x=ky =^ fix) =k f{y)- 

Using this notion, it is not hard to prove the following cancellation law up to fc-th order 
infinitesimals. 

Theorem 40. Let m G *M, n G N>o, j G N" \ {0} and a G M^g. Moreover let us consider 
/c G M defined by 

^ n 
1=1 

then 

1. \/h E -Dq^ X ■ ■ ■ X Da„ '■ ■ m = ■ Lkm 

2. If ■ m = for every h G Da^ x ■ ■ ■ x Da„, then m =k 
E.g. if n = 1 and ai = ji = 1 we have k = 2 and hence 

E D : h ■ m = h ■ L2Tn 

(yheD : h-m = 0) m =2 0. (68) 

Using (l68l) in (l64ll we obtain the final conclusion 

P-g^=2T--^ + G-p VxG (a,6) Vt G (0,+oo). (69) 



It is also interesting to note that not only small oscillations of the string implies (l69ll . 
but the converse is also true: the equation ( l69ll implies that necessary we must have small 
oscillations of the string, i.e. that tf{x,t) G -Dqo- Moreover, using the equality =2 up 
to second order infinitesimals, all the classical approximation tied with Hook's law, now 
become more clear. Indeed, we have the following 
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Theorem 41. Let a, b e R, with a < b; let 7 : [a,b] x [0, +00) — > 'M^, p : [a,b] x 
[0, +00) — > 'R and G, T : [a, b] x [0, +oc) — > 'R^ be non standard smooth functions and 
T E 'R be an invertible Fermat real. Let us suppose that the first component 71 of the curve 
is of the form 

7i(x,t) = [1 + ■x + /5(t) Vx,t, (70) 

with a{t) G l2- Then the unit tangent vector t{x,t) to the curve 7 exists and we can further 
suppose that the relations 

T{x,t) =2T ■t{x,t) (71) 
p-Sx- — ^ = T(x + Sx, t) - T(x, t) + G-p-5x, (72) 

holds for a every point (x, t) G (a, h) x [0, +00) and for every 6x G D. Finally, let us suppose 
that 

dip 

-;r-(x,t) is invertible. 
ox 

Then at this point (x, t) the following sentences are equivalent 

1. p(x, t) ■ ^(x, t) =2 T ■ ^(x, t) + G2(x, t) ■ p(x, t) 

2. (p{x,t) G I4. 

Finally, if ([2]) holds for every (x, t) G (a, 6) x [0,+oo), then 

length(7_,t) =36-0. 
To simplify the proof of this result, we need two lemmas. 

Lemma 42. Let a, b E R with a < b and let f , g : (a, b) — > 'R be non standard smooth 
functions such that 

fix) =2 g{x) Vx G (a, 6). 

Then 

/(x + /i) — /(x) = (7(x + h) — g{x) V/i G -D Vx G (a, 6) 

Lemma 43. Let m, h E *R, and suppose that m is invertible and < h < n, then the 
following properties are equivalent: 
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1. m ■ cos^ /i =2 m 

2. he h. 

Proof of Theorem [4Tt We firstly note that, assuming ( iTOl l. the tangent vector t{x,t) 
always exists in 'M. In fact we have ^{x,t) = 1 + a{t) so that both ^{x,t) and 
[^(''''^)]^ [^(''''^)]^ ^'■^ invertible; we can hence take its square root and then the 
inverse to define the unit tangent vector. Now we prove that ^ implies Let us take a 
generic 5x E D. Projecting (j72ll on 62 we get 



p ■ Sx 



T{x + 6x, t) -62- T(x, t) ■ 62 + G2 ■ p ■ Sx. 



But from (TfTl) and because smooth operations preserve =2, we get T ■ 62 =2 T ■ t ■ 62. 
Therefore, from Lemma [42] we obtain 



T(x + Sx, t) ■ 62 — T(x, t) ■ 62 = T • t(x + Sx, t) ■ 62 —T ■ t(x, t) ■ 62 

= T ■ sin Lp{x + 5x, t) —T ■ sin t) 



p ■ Sx 



T ■ sin (p{x + Sx, t) — T ■ sin (p{x, t) + G2 ■ p ■ Sx. 



(73) 



On the other hand, we can multiply ([T]) by Sx (so that =2 becomes =, see Theorem [40 
obtaining 



p ■ Sx 



T- 



d'~^2 d'~^2 
— {x + Sx,t)-—ix,t) 



+ G2 ■ p ■ Sx 



T ■ ta.nip(x + Sx,t) ■ -t—{x + Sx,t) — T tainp{x,t) ■ — — (x, t) + G2 ■ p ■ Sx, 

ox ox 



(74) 



Equating ([73l l and ( f74l l and canceling T we get 



simpix + Sx,t) — sin Lp(x,t) = tainpix + Sx,t) ■ -Tr—{x + Sx,t) — iaxiip{x,t) ■ -—{x,t) 

ox Ox 



dip 

ox ■ COS ip ■ — 

ox 



Sx 
Sx 
Sx 



1 


cos^ 




1 




cos^ 




1 





:^{x,t) + tany? ■ — -— 



d(p 



dip 



[l + a{t)] 



cos^ (y9 dx 



(75) 
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where, as usual, every function, if not otherwise indicated, is calculated at {x,t). Let us 
note that, in (TfSl) we have used the property 6x ■ a{t) = because 6x E D and a{t) G h; 
moreover, from ( IST! ) if (p = ^ we would have ^ ■ cosy9 = = ^ ■ siny? = 1 + a{t), which 
is impossible because a{t) G -Dqo- Setting, for simplicity, m := |^(x,t) G 'M, from ( 1751 ) and 
canceling 6x, we have 

m • cos^ (p =2 m, (76) 

By Lemma [43] this implies the conclusion. 

Vice versa, if ip is an infinitesimal of order less than or equal 4, then by Lemma [43] we 
obtain ([761 ) and we can go over again the previous passages in the opposite direction to prove 

©• 

Now, let us suppose that p{x,t) G I4 for every {x,t) G (a, 6) x [0, +cxd), then 



lenght(7_ 




dx 



1 2 



x,t) 



dx. 



(77) 



because a{t) G h and hence a{t)'^ = 0. But [1 + a{t)] ■ sin*/? = -^{x,t) ■ cosip, so 

^72 



9a: 



(x, t) = [1 + a{t)] tan 
= [l + a(t)] (p + 



p- 



= + y + «(^) ■ V^, 

because a;(t) G I2 and ip ^ I4 and hence Q;(t) ■ = 0. Substituting this in (TfTl) and using 



the derivation formula for the function a; t—>- vT+x we obtain 



l + 2a(t) + 



1 + - ■ <! 2a(t) + 



l + a{t) + - 



2 4 
1 + + Y + y + 



Therefore 
length(7_,i) 



dx 



h — a + a{t) ■ (6 — a) + 



'^{.x,tf ^ p{x,t) 



+ a{t) ■p{x,tf 



dx. (78) 
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Using the Theorem [38] it is not hard to prove that the last integral in (l78l) is an infinitesimal 
of order less than or equal 2, so the conclusion follows from the hypothesis a{t) G l2- 

m 

Proof of Lemma 142b First of all, from the hypothesis f{x) =2 g{x) for every x G {a,b), 
we get that 

°f{x) = °g{x) \fxe{a,b). (79) 
Now, let us fix a point x G (a, b). From Theorem [381 we obtain that we can write 

f(xi) = ao(xi) + ■ ai(xi) 

i 

g{xi) = bo{xi) + ^qj- bj{xi), 
j 

for every Xi G {x — 6,x + 6) C (a, 6) and where pi,qj G D^o and Oj, 6j are ordinary smooth 
functions defined in an open neighbourhood V of °x G (a, b) fl R. From (TfOl) we have 
ao(°a;i) = &o(°a^i) for every a^i G 'l^ so that = bo on V and hence also 'ao = on 'V^. 
Therefore 

/(r) - (7(r) = ■ a^(r) - J^g^ • &,(r) Vr G (a, 6) n R. (80) 

« i 

This difference must have order less than or equal 2 because /(r) =2 (?(?"), so 



UJ 



^Pi-ai{r) -^Qj -b 



maxuj [pi ■ ai{r)] V maxcj [g^ ■ 6j(r)] < 2. 



Let us suppose, for simplicity, that uj{pi ■ ai(r)) is this term of maximum order. Because 
ai(r) G M it must be that cij(pi) < 2 and hence a\souj{pi) <uj{pi) < 2 and co'(gj) <uj{pi) < 2. 
Finally we have 

f{x + h)-f{x) = h-f'ix) 

= h ■ ttQ^x) + h ■ Pi ■ a-(x), 

i 

but aQ{x) = &o(x) because ao = &o and h ■ Pi = because uj{h) < 2 and uj{pi) < 2; we hence 
obtain 

fix + h)-f{x) = h-b'oix) 

= h-b'o{x) + ^h-qj-b'j{x) 

= h ■ g\x) 

= g{x + h)- g{x). 
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Proof of Lemma I43t If m ■ cos^ h =2 m, then the standard parts of both sides must be 
equal 



(m ■ cos^ h) 



m ■ cos 



3 /o 



m 



m. 



By hypotheses m is invertible, hence °m 7^ and we obtain that °h = because < h < n, 
i.e. h G -Doo- Moreover, from infinitesimal Taylor's formula applied to cos/i, and from 
m ■ cos^ h =2 m we obtain 



m ■ 



1- E ' 



h 



2i 



\ 



--2 m 



m ■ (1 + a^h^ + 3a/i^ + 3a^/i^) =2 m 



m 



(1 + a ■ /i^) =2 m 



(2*)! 



G 'M and a := 3a^ + 3a + a^/i^ are invertible Fermat 



where a := — 

reals. From this we get m-a-Zi^ =2 and hence /i^ =2 0, i.e. c<j(/i^) < 2 and hence ujiK) < 4. 

Vice versa, if h is an infinitesimal of order less than or equal 4 (so that (p"' = if n > 5) 
we have 



cos h 



2^4! 

3— + 3—. 
2 4! 



Therefore, m ■ cos^ h = m — 3mh'^ ■ ^| — 3|^ j so that m ■ cos^ h 
is an infinitesimal of order a;(/i^) < 2, i.e. mcos^ /i 



-3m/i^ 



=2 



The reader with a certain knowledge of SDG had surely noted that this deduction of the 
wave equation cannot be reproduced in SDG because of the use of non standard smooth 
functions, of the use of equalities up to k-th order infinitesimals and because of the frequent 
use of the useful Theorem [12] to decide products of powers of nilpotent infinitesimals. 
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XV. CONCLUSIONS 

The problem to turn informal infinitesimal methods into a rigorous theory has been faced 
by several authors. The most used theories, i.e. NSA and SDG, require a good knowledge 
of Mathematical Logic and a strong formal control. Some others, like Weil functors (see 



e.g. Kriegl and Michor [2l|) or the Levi-Civita field (see e.g. Shamseddine [25|) are mainly 
based on formal/algebraic methods and sometimes lack the intuitive meaning. In this initial 
work, we have shown that it is possible to bypass the inconsistency of SIA with classical 
logic modifying the Kock-Lawvere axiom (see e.g. Lavendhomme [2^) and keeping always a 
very good intuitive meaning. We have seen how to define the algebraic operations between 
this type of nilpotent infinitesimals, infinitesimal Taylor formula and order properties. In 
the final part we have seen several elementary examples of the use of these infinitesimals, 
some of them taken from classical deductions of elementary Physics. In our opinion, these 
examples are able to show that some results that frequently may appear as unnatural in a 
standard context, using Fermat reals can be discovered, even by suitably designed algorithm. 
Moreover, our generalization of the classical proof of the wave equation have shown that a 
rigorous theory of infinitesimals permits to obtain results that are not accessible using only 
an intuitive approach. 
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